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Abstract 


A  system  to  measure  the  non-steady  velocities  in  the  model  rotor 
wake  was  modified  to  provide  velocity  data  in  a  form  more  useful  to  the 
user  of  the  data.  Hot  wire  angle  parameters  were  modified  so  as  to 
provide  better  correlations  by  using  cosine  law  in  place  of  the  usual 
K  technique  common  in  hot-wire  anemometry.  The  results  of  the  non¬ 
steady  measurements  were  presented  in  the  form  of  motion  pictures  of  the 
non-steady  computer  generated  vector  plots. 
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CHAPTER  I 


INTRODUCTION 

Helicopter  rotor  wake  characteristics  during  low  speed  terrain 
flight  have  become  the  subject  of  increased  interest  in  recent  years, 
particularly  in  military  applications. 

The  interaction  of  blade  tip  vortices  at  transition  speeds, 

(0  <  u  <  0.2)  and  the  resulting  wake  patterns  are  significant  in  terms 
of  their  effect  on  aerodynamic  surfaces  such  as  wings,  stabilizers,  or 
fuselage.  The  apparent  inability  of  theoretical  methods  to  accurately 
predict  wake  velocities  in  the  transition  flight  regime  has  led  to 
many  experimental  investigations.  Refs.  (5,6,8). 

This  work  is  an  extension  of  the  experimental  investigation  of 
Ref.  (8),  in  which  the  measurement  capability  of  non-steady  rotor  wake 
velocities  was  established.  Presented  here  are  the  computer  based 
non-steady  rotor  wake  velocity  data  reduction  and  presentation 
techniques. 

The  massive  amounts  of  data  being  manipulated  necessitated  the  use 
of  a  computer  facility.  The  system  used  here  consists  of  a  DEC  PDP 
11/60  minicomputer  and  a  GENRAD  analog  to  digital  converter.  Flow 
measurements  were  made  in  a  4  x  8  ft.  wind  tunnel,  of  a  2-bladed  model 
helicopter  rotor  using  a  3-dimensional  hot  film  anemometer  system. 
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CHAPTER  II 


EXPERIMENTAL  APPARATUS 


2.1  The  Test  Rig 


The  data  for  this  work  was  acquired  using  the  Mechanical  Engi¬ 
neering  Department  4  ft.  x  8  ft.  wind  tunnel  facility  at  The  Ohio 
State  University.  The  test  section  contained  a  2.5  ft.  dia.,  two- 
bladed,  teetering  rotor.  Collective  pitch  and  rotor  shaft  tilt  angle 
were  adjustable  but  once  set  for  a  test,  they  remain  constant.  The 
rotor  details  are  listed  in  Table  2.1. 

Table  2.1.  Details  of  Rotor 

Number  of  Blades  2 

Rotor  Radius,  ft.  1.25 

Blade  Chord,  in.  2.1875 

Rotor  Solidity  bc/irR  .0928 

Root  Cutout  »R  12.08 

81ade  Taper  Ratio  1 

Blade  Twist  8°  (from  root  to  tip) 

Airfoil  Section  NACA  0012 

Rotor  speed  is  variable  up  to  a  maximum  of  3000  RPM  yielding  a  maximum 
tip  speed  of  400  ft/sec.  The  wind  tunnel  is  capable  of  stream  veloci¬ 
ties  up  to  a  maximum  of  75  ft/sec.  Figure  2-1  shows  the  test  section 
and  the  rotor  orientation  within  the  test  section. 
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Figure  2-1.  Rotor  Location  in  Test  Section  Looking  Downstream 

2.2  The  Hot-Film  Anemometry 

The  rotor  induced  velocities  were  measured  by  means  of  a  constant 
temperature,  3-dimensional  hot-film  anemometer  system  incorporating  a 
TSI  model  1294-60-18  probe,  and  3  DISA  55005/1 02C  control  circuits. 

The  DC  power  supply  for  the  control  circuits  and  signal  amplifiers  is  a 
SOLA  +  15V,  1.6A  model  83-15-2216. 

This  anemometry  system  used  in  these  tests  differs  from  that  of 
Ref.  (8  )  in  that  changes  became  necessary  in  order  to  eliminate  severe 
calibration  drift  problems  encountered  at  the  onset  of  this  work.  The 
system  as  used  here  is  shown  in  Figure  2-2.  The  differences  between 
this  and  the  system  of  Ref.  (8  )  include:  1)  the  elimination  of  lead 
compensators;  2)  the  use  of  .006  in.  diameter  probe  as  opposed  to  .002 
in.  diameter;  3)  the  substitution  of  3  two  stage  integrated  circuit 
signal  amplifiers  for  the  3  DANA  2210  amplifiers  and;  4)  minor  modifi¬ 
cation  of  the  DISA  control  circuits.  The  lead  compensators  were 
removed  for  two  reasons.  First,  it  was  suspected  that  they  may  have 
been  responsible  for  the  drift  phenomenon.  Though  they  may  have 


uibAaMiriil 


j 


TSI  1294-60-18 


Figure  2-2.  Block  Diagram  of  the  Anemometer  System 
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contributed  to  this  problem,  their  removal  did  not  solve  the  problem. 
Secondly,  their  effectiveness  in  extending  frequency  response  was  being 
negated  by  frequency  limiting  software  techniques  which  were  employed 
in  an  effort  to  reproduce  results  previously  obtained.  Details  of  this 
software  approach  are  included  in  Appendix  A. 

The  three  integrated  circuit  amplifiers  provided  the  variability 
needed  in  controlling  the  output  voltage  range  for  the  range  of  input 
velocities,  0  to  65  ft/sec.  The  range  of  amplification  is  continuously 
variable  from  a  factor  of  1.0  to  100.  The  DANA  units  have  discreet 
amplification  factors  making  it  necessary  to  add  voltage  dividers  for 
"fine  tuning".  This  is  particularly  cumbersome  since  changing  probes 
necessitates  a  change  in  the  voltage  divider.  Figure  2-3  shows  the  new 
amplifier  circuit.  Note  here  also  that  a  bias  voltage  adjustment 
potentiometer  is  provided. 

After  an  attempt  at  designing  a  completely  new  anemometer  control 

circuit  met  with  only  moderate  success,  see  Appendix  B,  the  existing 

DISA  control  circuitry  was  repaired.  All  transistors  in  Channel  A 

were  replaced,  see  Figure  2-4,  as  was  the  TD101  transistor  pair  in 

Channel  C.  It  was  also  discovered  that  the  capacitor  between  the  base 

of  the  2N3906  transistor  and  ground  has  a  significant  effect  on  the 

frequency  response  and  output  drift  of  the  system.  After  trying 

several  values  in  the  range  between  ,056yF  and  .22pF  a  value  of  .15pF 

was  chosen.  The  other  modification  was  the  addition  of  a  100Q  OHMITE 

model  4092  resistor  in  series  with  the  IK  potentiometer  in  the  bridge. 

This  was  necessary  to  accommodate  a  probe  which  had  an  operating 

resistance,  R.„,  equal  to  10.1 3C2.  Since  the  overheat  ratio  is  obtained 

op 


l*Z9ZA 
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by  multiplying  the  sum  of  Rgp  and  R^eacj  wire  by  100,  the  resistance  in 
this  leg  of  the  bridge  must  be 


RP0T 


rpot 


=  (R0P  +  Rlw)  x  100 

=  (10.13  +  .30)  x  100 
=  1043fi 


With  the  addition  of  the  lOOfi  series  resistor  this  value  of  RpQj  can  be 
reached. 

The  anemometer  circuitry  as  described  here  should  be  considered 
only  marginally  adequate  for  the  task  being  performed.  A  new 
circuit  of  the  same  design  would  improve  the  reliability  and  stability. 


2.3  Data  Acquisition 

2.3.1  The  Probe  Traverse  Mechanism 

The  probe  is  mounted  on  a  traverse  mechanism  which  allows 
the  operator  to  position  the  probe  at  any  point  in  the  test  volume. 
The  spatial  orientation  of  the  probe  in  the  tunnel  is  depicted  in 
Figure  2-5.  The  traverse  mechanism  employs  motorized  positioning  in 
the  Y  and  Z-directions  and  a  manually  operated  pully  arrangement  for 
the  X-direction.  The  test  volume  measures  68  in.  x  40  in.  x  15  in., 
(X,  Y,  and  Z  respectively).  In  the  X-direction  data  is  taken  at  each 
4  in.  increment  starting  at  X  =  28  in.  going  to  X  =  -  40  in.  In  the 
Y-direction  data  is  taken  at  4  in.  increments  starting  at  Y  =  20  in. 
going  to  Y  *  -  20  in.  In  the  Z-direction  data  is  taken  at  3  in. 
increments  starting  at  Z  *  12  in.  going  to  Z  *  -  3  in.  Thus  the  test 
volume  for  a  complete  data  set  contains  a  matrix  of  1188  points. 
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Figure  2-5.  The  Spatial  Orientation  of  the  Probe  in  the  Wind  Tunnel 
However,  due  to  rotor  interference,  147  points  cannot  be  measured  as 
indicated  by  the  shaded  region  in  Figure  2-6.  Note  also  that  the 
rotor  hub  center  is  the  point  <0,0, 0>  for  a  shaft  tilt  angle  equal  to 
zero. 

2.3.2  The  Computer  Interface 

The  anemometer  system  is  linked  to  the  DEC  POP  11/60  mini¬ 
computer  via  an  interface  circuit  and  a  GENRAD  A/D  converter.  The 

interconnection  of  the  two  systems  is  shown  in  Figure  2-7.  The  inclu¬ 
sion  of  the  delay  lines  in  Channels  B  and  C  is  to  allow  for  the  read 

time  cycle  of  the  A/D.  Without  the  delay  lines.  A,  B,  and  C  would  be 
read  at  slightly  different  points  in  the  rotor  revolution. 

The  interface  circuit  is  actually  a  pulse  shaping  circuit  which 
triggers  the  A/D.  The  Interface  receives  pulses  from  two  magnetic 
pickups  mounted  on  the  rotor  support  structure.  One  pickup  senses  each 


tooth  of  a  72  tooth  gear  turning  with  the  rotor.  The  second  pickup  is 


Figure  2-6.  Rotor  Interference  Volume  Relative  to  Test  Volume 
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a  once  per  revolution  signal.  The  output  of  the  interface  is  a  pulse 
train  which  triggers  the  A/D  to  sample  sensors  A,  B,  and  C  72  times 
per  revolution  (every  5  degrees)  for  as  many  revolutions  as  required. 

In  this  case  14  revolutions  were  used  giving  a  total  of  72  x  14  x  3 
or  3024  pieces  of  data  at  each  point  in  the  test  volume.  Assuming  then 
that  the  signal  is  periodic,  an  average  is  obtained  so  that  only  72  x  3 
or  216  pieces  of  data  are  actually  stored  in  computer  memory.  A 
further  reduction  in  the  amount  of  memory  space  needed  is  accomplished 
through  a  double  precision  integer  coding  scheme.  Essentially  this 
technique  codes  the  velocity  Vx  in  the  first  3  digits  of  a  9-digit 
integer,  Vy  in  the  second  3  and  V2  in  the  last  3.  This  uses  less 
memory  space  than  3  separate  real  numbers. 


CHAPTER  III 


PROBE  DIRECTIONAL  SENSITIVITY  AND  CALIBRATION 


3.1  Directional  Sensitivity  of  Cylindrical  Sensors 

The  output  of  the  hot-film  anemometer  is  a  voltage  which  is  relat¬ 
ed  to  the  fluid  velocity  as: 

E2  =  A  +  BU1/2  (3-1) 

where  U  is  the  stream  velocity,  E  is  the  output  voltage,  and  A  and  B 
are  constants  determined  by  the  circuitry  and  the  physical  properties 
of  the  sensor  and  fluid  medium.  However,  if  the  fluid  velocity  is  not 
perpendicular  to  the  sensor,  equation  3-1  does  not  hold.  To  overcome 
this  problem,  an  "effective  cooling  velocity"  can  be  defined  as  the 
perpendicular  velocity  which  produces  the  same  cooling  effect  as  the 
actual,  non-perpendicular  velocity.  This  effective  cooling  velocity 
must  then  be  correlated  to  the  true  stream  velocity.  The  relationship 
used  in  Ref.  (6)  is  given  as: 

V2ff  *  U2(cos2$  +  K2sin2o)  (3-2) 

where: 

Vgff.  =  effective  cooling  velocity 

U  =  stream  velocity 

$  *  angle  between  the  free  stream  vector  and  the  normal  to 
the  sensor.  See  Figure  3-1. 

K  *  dimensionless  constant 
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2  2 

The  term  K  sin  <p  accounts  for  the  cooling  effect  of  flow  along  the 
sensor. 

Another  correlation  which  is  used  is  that  of  Ref.  given  by: 

Veff  =  U  cosP  ♦  (3-3) 

where  the  exponent,  P,  is  used  to  account  for  cooling  effects  of  flow 
along  the  sensor.  This  relationship  is  used  here  in  order  to  simplify 
the  data  reduction  procedure. 

3.2  Method  of  Data  Reduction 

The  aim  of  the  data  reduction  procedure  is  to  derive  from  the 
voltage  outputs  of  each  sensor  the  spatial  orientation  of  the  stream 
vector  in  the  tunnel.  It  is  therefore  critically  important  that  the 
angles,  a,  6,  and  y  of  Figure  3-2  be  accurately  determined. 

Using  the  method  of  3.1,  each  of  the  3  mutually  perpendicular 
sensors  will  have  an  equation  for  the  effective  cooling  velocity  across 
it. 

VfA  *  VA  '  U  cosP  «A  <3-4> 

VeffB  ■  VB  ■  U  c°sP  *8  (3-5) 

VeffC  "  vc  *  u  cos  t3'6) 

Since  V^,  Vg,  V^,  and  P  (note  that  p  may  be  different  for  each  sensor, 
but  are  used  as  the  same  value)  will  be  known  from  the  calibration,  U 
must  now  be  found  in  terms  of  these  quantities. 

In  order  to  make  use  of  the  geometric  constraint  that: 
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COS2<J>a  +  COS2<t>g  +  COS2$c  =  2 
equations  3-4  to  3-6  are  arranged  in  the  form 
V.  1/P 

(-jj-)  =  cos-fK  i  =  A,B,C 

It  then  follows  that 


v  i/p  e 
*[(-£)  ]  »  2 

and  in  expanded  form  this  becomes 

h  ^  -  z 


u 


which  can  be  solved  for  U 


U  =  [- 


2/P  +  y  2/P  +  2/P  P/2 


-] 


(3-7) 


(3-8) 


(3-9) 


(3-10) 


(3-11) 


Now  in  order  to  find  the  angles  <J>A,  d»B ,  and  ^  which  will  locate  U 
in  the  probe  system  of  coordinates,  equations  3-4  to  3-6  will  be  solved 
for  these  angles. 


*  cos'1 

V  1/P 

A 

1  u ' 

(3-12) 

4>g  =  cos'1 

Vg1/P 

<<T> 

(3-13) 

-1 

a  cos 

V  l/p 

(if) 

(3-14) 

With  4>g  and  known,  the  vector  U  has  been  located  in  probe  coordi¬ 
nate  system. 


3.3  Transformation  of  Velocities  from  the  Probe  to  the  Tunnel  Coordi 


nate  System 

Due  to  the  nature  of  the  probe  geometry  and  its  orientation  within 
the  tunnel,  the  vector  U  must  be  transformed  into  its  X,  Y,  and  Z 
components  in  the  tunnel  coordinate  system  shown  in  Figure  3-3. 

Though  the  probe  axis  is  parallel  to  the  tunnel  Z  direction  the  indi¬ 
vidual  sensors  are  each  inclined  at  an  angle  of  54.74°  to  the  probe 
axis.  Figure  3-3  also  shows  the  vectors  OA,  OB,  and  OC  which  repre¬ 
sent  flow  in  the  direction  of  sensors  A,  B,  and  C  respectively.  Note 
also  that  the  vector  OC  lies  in  the  X-Z  plane. 

If  UA,  Ug,  and  Uc  are  defined  as  the  components  of  U  in  the 
direction  of  sensors  A,  B,  and  C  respectively  their  equations  are: 


UA  = 

Usin$A 

(3-15) 

UB 

Usintfg 

(3-16) 

UC  = 

Usimj>£ 

(3-17) 

where  it  should  be  recalled  that  is  defined  as  the  angle  between  U 
and  the  normal  to  the  sensor. 

From  the  geometry  of  the  probe  and  its  location  with  respect  to 
X,  Y,  and  Z  in  Fig.  3-3  the  angles  EOC,  EOD,  AOD  and  BOD  are: 

<E0C  =  35.26° 

<E0D  =  54.74° 

<A0D  =  45° 

<B0D  =  45° 

Using  these  angles  the  transformation  equations  3-18  through  3-20  are 
defined  as 
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Vx  =  -  0.40825(Ua  +  UB)  -  0.8165  Uc  (3-18) 

Vy  =  0.7071  (UA  -  Ug)  (3-19) 

Vz  =  0.57735(Ua  +  Ug  -  Uc)  (3-20) 

The  quantities  V  ,  V  ,  and  V  are  the  components  of  the  vector  U  in  the 
x  y  z 

tunnel  coordinate  directions  x,  y,  and  z. 


3.4  Probe  Calibration 


3.4.1  Velocity  Calibration 

The  velocity  calibration  of  the  probe  was  performed  using 
a  calibration  air  chamber  with  an  exit  nozzle  diameter  of  0.75  in.  The 
chamber  static  pressure  is  measured  using  an  inclined  manometer  from 
which  the  nozzle  exit  air  velocity  is  determined.  Each  sensor  is  in 
turn  positioned  such  that  its  support  prongs  lie  in  the  horizontal 
plane  and  the  sensor  is  within  0.5  in.  of  the  vertical  exit  plane  of 
the  nozzle  and  parallel  to  it.  The  jet  velocity  is  subsequently  varied 
from  0  to  65  ft/sec.  as  the  output  voltage  is  recorded  for  each  veloci¬ 
ty.  A  total  of  13  velocities  are  used  in  the  range  from  0  to  65  ft/sec 
and  a  canned  computer  program  is  used  to  fit  a  fourth  order  polynomial 
through  these  points.  The  velocity  calibration  curve  for  channel  A 
used  in  this  work  is  shown  as  an  example  in  Fig.  3-4.  The  equations 
for  this  curve  and  similar  curves  for  B  and  C  are  entered  in  the  com¬ 
puter  program,  C0SPX.FTM  (Appendix  C),  and  are  of  the  form 


V  =  C  E 
veff  V 


+  c3e3  +  c2e2  +  ^E  +  Cq 


(3-21) 


3.4.2  Directional  Sensitivity  Calibration 

In  previous  work  the  directional  sensitivity  of  the  probe 

2 

sensors  has  been  accounted  for  by  use  of  the  constant,  K  ,  in  equation 
3-2.  This  work  employs  the  constant,  p,  of  equation  3-3  as  a  result  of 
calibration  data  which  suggested  that  this  method  would  work  well. 

The  angular  calibration  data  taken  is  independent  of  which  method 
is  used.  The  probe  is  positioned  on  the  calibrator  in  the  same  manner 
as  in  the  velocity  calibration,  only  now  the  angle  <p  (see  Fig.  3-1) 
is  varied  from  0°  to  65°  in  5°  increments.  The  angle  $  was  not  varied 
beyond  65°  due  to  suspected  interference  of  the  sensor  support  prongs. 
The  output  voltage  at  each  angle  is  recorded  and  this  process  is 
repeated  for  8  velocities  in  the  range  from  10  ft/sec  to  65  ft/sec. 

At  this  stage  a  computer  program  is  used  to  calculate  Veff/U 
for  each  angle  <j>  using  equation  3-21.  It  was  decided  to  fit  the 
equation. 


eff  . 
U 


cos 


(3-22) 


to  the  V  ,,/U  data  at  the  4>  =  65°  point  by  choosinq  an  appropriate 
eTT  v 

veff 

value  for  p.  The  calibration  data  indicated  that  — g—  at  <j>  =  65°  was 
nearly  a  constant  for  all  values  of  U.  This  being  the  case,  a  constant 
value  of  p  can  be  used.  The  values  of  p  calculated  for  each  sensor  are 


Sensor  A 
Sensor  B 
Sensor  C 


p  =  0.76 

p  =  0.80 

p  =  0.78 


Note  here  that  an  average  value  of  p  must  be  used  in  equation  3-11. 

Figure  3-5  shows  the  correlation  between  the  calibration  data  (for 
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sensor  C  at  li  =  16  ft/sec  and  U  =  32  ft/sec)  and  the  cosP<{>  and 

2 

methods  with  p  =  0.78  and  K  =  0.10.  It  was  observed  from  this  plot 
that  U  has  little  effect  on  the  variation  of  Vg^/U  with  $.  This  is 
most  likely  due  to  the  relatively  large  diameter  sensors  in  use  here. 


Figure  3-3. 


Probe  Coordinate  System  Relative  to  Tunnel  Coordinate 
System 


Correlation  Between  Calibration 
Methods 


CHAPTER  IV 


PRESENTATION  AND  DISCUSSION  OF  TEST  RESULTS 

4.1  Test  Conditions 

The  anemometer  system  of  section  2.2  was  used  to  acquire  the 
averaged  instantaneous  data  for  an  entire  case  which  consists  of  1138 
points  of  velocity  information.  The  frequency  response  of  the  system 
as  determined  by  the  square  wave  test.  Ref.  (3  ),  is  approximately 
400  Hz. 

The  parameters  of  the  flight  condition  were  as  follows: 


Blade  Tip  Speed,  ft/sec  300 
Advance  Ratio  0.06 
Collective  Pitch  (75^R),  degrees  8 
Blade  Twist,  degrees  8 
Rotor  Shaft  Tilt  Angle,  degrees  8 


4.2  Data  Presentation 

The  velocity  information  is  presented  here  in  three  forms. 

Figures  4-1  through  4-3  show  the  y-z,  x-z  and  x-y  vectors  respectively 
for  the  blade  azimuth,  ip  s  355°.  Figures  4-4  through  4-6  show  a 
selected  x,  y  and  z  plane  respectively,  for  blade  angles  from  0°  to 
315°  in  45°  Increments.  The  data  is  available  in  5°  increments  but  is 
not  presented  as  such  in  the  interest  of  brevity.  It  should  be  noted 
here  that  In  Figures  4-1  to  4-6  the  ordinate  and  abscissa  scales 
represent  the  non'-dimensionalized  tunnel  location  and  the  "+"  signs 
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represent  the  origin  of  the  vectors.  Velocity  magnitude  scaling  can  be 

done  using  the  scale  provided  in  the  legend  of  each  plot. 

The  third  format  used  to  present  the  data  is  that  of  Figure  4-7. 

These  graphs  show  the  variation  of  velocity  (V  ,  V  and  V  )  in  one 

x  y  z 

rotor  revolution  for  an  individual  point  in  the  test  volume.  A 
traverse  of  six  points  in  the  z-direction  is  depicted  in  Figure  4-7, 
starting  at  Z/R  =  0.8  and  ending  at  Z/R  =  -0.2.  Table  4-1  gives  the 
non-dimensional ized  tunnel  location  and  point  number,  where  the  point 
number  is  in  reference  to  the  1138  point  matrix  of  the  test  volume. 
Table  4-1.  Tunnel  Location  Reference  for  Fig.  4-7 


POINT  # 

X/R 

Y/R 

Z/R 

444 

-1.07 

0.27 

0.8 

462 

-1.07 

0.27 

0.6 

480 

-1.07 

0.27 

0.4 

498 

-1.07 

0.27 

0.2 

516 

-1.07 

0.27 

0.0 

534 

-1.07 

0.27 

-0.2 

The  vector  plots  of  Figures  4-1  through  4-6  were  made  using  the 
program  VECTOR. FTN  whereas  the  plots  of  Figure  4-7  were  made  using  the 
program  EZPLOT. FTN  in  conjunction  with  the  system  graphing  routine 
called  EASY  GRAPHING,  see  Appendix  C.  Also  in  Appendix  C  are  the 
tabular  data  for  the  plots  of  Figure  4-7. 

The  data  of  the  azimuth- varying  vector  plots  of  Figures  4-4 
through  4-6  are  also  available  in  an  animated  film  version.  The  film¬ 
ing  technique  is  presented  in  Appendix  D. 
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Figure  4-1.  continued 
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Figure  4-1.  continued 
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Figure  4-3.  continued 


-:.4 


—  20  FT/S3C 


X/S  -- 1.07 


Figure  4-4. 


ip=0° 


'jJ-45a 


ip-30 


'i-l  Vectors  for  X/R  *  -1.07  at  Several  Blade  Angles 


Figure  4-4.  continued 


Figure  4-4.  concluded 


MSN  l>\N  WSN  HSH 


Figure  4-5. 


0.4 

o.o 

0.4 

0.8 


0.4 

0.0 

0.4 

0.8 


0.4 

0.0 

0.4 

0.8 


0.4 

0.0 

0.4 

0.8 


— 20  FT/SIC  Y/R  -  0.0 

*  - — . k'\\\V\ 

. \UU\\ 

- -  \  \  \ 

v 

- 

■  ■ _ i-  —i _ i  i _ 

Y/R  ■  0.0 

*  - . \\\\\^ 

. \u\\\\ 

"" . \\\\\\\ 

- - -  -  -  s  *-  \  \  V 

v  K  s  * - 'v’sV^V.  v  v 

i _ i _ i _ i _ 

Y/R  -  0.0 

’•  *  •  •  ♦  ♦  ♦  *•  v  \  \  \  \  S, 

— .  v\\\\\v 

. \\\  \\\\ 

*"  -  ^  ~  “  -  *'\s-\\\\  \  \  \ 
r'~‘ - - 

V. V-sS^V 

i _ I _ I _ I - 1 - 1 - 

20  PT/SSC  Y/R  -0.0 

. . fc\\\\\^ 

. \u\\\\ 

\\U\\N 

\\W  \  \  \ 

'  ^VSjVO>V 

■  1  1 

1.6  0.8  0.0  -0.8  -1.6  -2.4 

FORWARD  X/*  AFT 


\p-0° 


\p-4S 


\p=90 


\p-l3i 


X-Z  Vectors  for  Y/R  ■  Q.O  at  Several  Blade  Angles 


m.+i 


Figure  4-7.  Azimuth  Varying  Velocities  for  1  Blade  Revolution 
Z  Traverse  at  Y/R  *  0.27  and  X/R  *  -1.07 
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4.3  Discussion  of  Results 

As  can  be  seen  by  the  representative  data  in  Figures  4-1  through 
4-7,  there  is  a  vast  amount  of  information  available  for  analysis. 
Ideally,  the  data  should  pertain  to  a  helicopter  in  level  unaccelerated 
forward  flight.  This  condition  exists  when  the  thrust  vector  is  coin¬ 
cident  with  the  rotor  shaft  axis.  However,  no  provision  was  made  to 
ensure  this  condition.  Consequently,  no  specific  flight  attitude  can 
be  presumed.  Despite  this  limitation,  the  data  presented  demonstrates 
the  feasibility  of  obtaining  qualitative  and  quantitative  averaged 
instantaneous  rotor  wake  velocity  information. 

Of  particular  interest  in  the  transition  flight  regime  is  the 
nature  of  the  rollup  phenomenon  of  the  trailing  vortices.  This  is 
shown  clearly  in  Figure  4-1  beginning  dt  X/R  *  0.27  suggesting  that 
the  vorticity  sheet  rolls  up  very  quickly.  Note  also  that  this  activi¬ 
ty  is  occurring  first  on  the  advancing  side.  At  location  X/R  *  -1.07 
the  rollup  is  fully  developed  and  a  very  strong  downwash  velocity  is 
exhibited.  Once  again  the  plots  fail  to  show  axisymmetric  flow 
patterns  in  that  the  downwash  velocities  on  the  advancing  side  are 
greater  in  magnitude  than  on  the  retreating  side.  Nonuniform  induced 
velocities  suggest  a  nonuniform  disk  loading.  Nonuniform  disk  loading 
is  expected  for  a  trimmed  rotor,  but  since  the  rotor  has  not  been 
trimmed  it  is  impossible  to  attribute  the  nonuniformity  to  a  specific 
cause.  The  anticipated  downward  and  aft  motion  of  the  trailing  vortex 
center  can  be  seen  when  comparing  the  apparent  center  on  the  advancing 
side  for  X/R  *  -1.07  and  X/R  =  -2.67.  This  downward  progression  is  not 
well  defined  on  the  retreating  side. 
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The  azimuth- varying  Y-Z  vector  plots  of  Figure  4-4  do  not  show 
obvious  changes  in  all  vectors  from  one  blade  position  to  the  next, 
however  upon  closer  inspection,  some  specific  vectors  are  widely  vary¬ 
ing.  Take  as  an  example  the  vector  at  the  point  Y/R  =  0.53  and 
Z/R  =  0.0.  During  the  blade  revolution  its  magnitude  changes  by  a 
factor  of  4  and  its  direction  changes  +  20°  from  the  vertical  direc¬ 
tion. 

The  lateral  progression  of  X-Z  planes  (Fig.  4-2)  from  the 
retreating  to  advancing  sides  shows  the  expected  strong  downwash 
patterns  typical  of  flows  for  low  advance  ratios.  The  downwash  is  most 
pronounced  at  the  Y/R  =  +0.53  locations.  At  the  location  Y/R  =0.0 
the  effect  of  the  tip  vortices  rotating  in  the  X-Z  plane  (this  will  be 
the  case  at  i|;  =  0°)  is  reflected  in  the  behavior  of  the  X-component. 

The  counterclockwise  rotation  of  a  vortex  in  the  region  just  below  and 
just  aft  of  the  rotor  will  add  to  the  stream  velocity  in  the  vicinity 
of  Z/R  =  0.4  to  0.8  and  subtract  from  it  in  the  0.0  to  -0.2  region. 
Beyond  the  blade  tip  path  in  both  directions  along  the  lateral  axis 
and  just  forward  of  the  rotor  these  same  tip  vortices  cause  a  slight 
upflow  which  can  be  seen  clearly  in  Fig.  4-2  at  locations  Y/R  =  +  1.33 
and  Y/R  =  0.0. 

A  very  obvious  indication  of  the  azimuth-wise  variation  of  the 
flow  can  be  seen  in  Fig.  4-5.  Bear  in  mind  that  at  ip  *  0°  the  blades 
are  in  the  plane  of  the  page  and  rotating  clockwise  as  viewed  from 
above  the  rotor  disk.  Noting  the  two  vectors  with  circled  origins  in 
the  <|«  »  0°  plot  and  tracking  them  through  the  revolution  it  can  be  seen 
that  the  X-component  of  these  two  vectors  in  the  ip  =  45°  plot  now  show 


the  influence  of  a  clockwise  rotating  tip  vortex.  At  <p  =  135°  the 
X-components  are  nearly  equal  again,  followed  by  a  nearly  identical 
repetition  of  the  vortex  influence  at  p  =  225°. 

The  X-Y  plots  of  Fig.  4-3  show  a  near  total  negation  of  the  stream 
velocity  at  Y/R  =  0.0,  Z/R  =  -0.2  and  X/R  =  -1.07  to  -2.4,  again  show¬ 
ing  the  tip  vortex  flow. 

The  azimuth- varying  velocities,  VX,  VY  and  VZ  of  Fig.  4-7  substan¬ 
tiates  what  has  been  shown  in  Figs.  4-1  to  4-6.  Note  the  clear  two  per 
revolution  frequency  appearing  in  the  X  and  Z  velocities  for  points 
516  and  534.  Moving  further  below  the  rotor  to  points  498,  480,  462 
and  444  shows  a  degeneration  of  this  frequency  probably  caused  by 
the  probe  being  more  removed  from  intense  vortex  activity. 


CHAPTER  V 


CONCL.  ONS  AND  RECOMMENDATIONS 

In  its  present  form,  the  entire  system,  including  the  computer 
programs,  is  an  invaluable  tool  in  the  data  acquisition  and  presenta¬ 
tion  process.  The  use  of  the  Mechanical  Engineering  Department  mini¬ 
computer  and  A/D  converter  greatly  expanded  the  capability  and  effi¬ 
ciency  of  the  system.  The  time  needed  to  take  the  data  for  the  case 
was  approximately  8  hours.  This  is  an  improvement  over  the  previous 
method  which  took  20  hours. 

The  calibration  process  is  currently  very  cumbersome  and  time 
consuming.  It  is  suggested  that  the  A/D  converter  and  the  minicomputer 
be  used  instead  of  the  current  method  of  using  a  digital  voltmeter  and 
the  IBM  470  batch  system.  This  would  eliminate  hours  of  work  with 
computer  cards.  The  calibration  results  could  be  stored  in  a  direct 
access  file  which  would  subsequently  be  read  when  the  data  taking 
program  is  run  (i.e.  COSPX.FTN).  This  would  eliminate  the  need  to 
enter  the  constants  of  the  polynomials  by  hand. 

Also  of  interest  in  terms  of  streamlining  the  tunnel  operators 
job  would  be  to  link  the  traverse  mechanism  to  the  minicomputer  which 
would  automatically  control  the  position  of  the  probe  in  the  tunnel. 

A  total  redesign  of  the  traverse  mechanism  and  the  purchase  of  a 
microprocessor  would  be  necessary.  The  microprocessor  would  allow  the 
flexibility  of  programmability  in  the  event  the  test  volume  should  be 
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changed.  In  contrast,  a  hard  wired  circuit  to  perform  the  same  task 
would  require  extensive  rework  if  any  test  volume  changes  are  made. 
The  new  traverse  mechanism  design  should  allow  complete  motion  from 
the  floor  to  the  ceiling  of  the  tunnel  without  the  currently  required 
operator  assistance. 
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CONCLUSIONS 

The  filmed  animation  of  the  rotor  flowfield  for  a  single  rotor 
revolution  demonstrates  the  feasibility  of  this  format  for  flow 
visualization.  The  overall  trends  of  the  data  show  agreement  with 
time-averaged  data  and  smoke  visualization  data  of  the  literature. 

Ref.  (5).  Expected  phenomenon  such  as  the  rollup  of  the  trailing 
vortex  in  the  Y-Z  planes  and  the  strong  downwash  velocities  just  aft  of 
the  rotor  are  both  typical  of  the  helicopter  rotor  in  forward  flight. 

Although  the  instantaneous  nature  of  the  measured  velocities  in 
the  x,  y  and  z  directions  is  readily  apparent  in  the  azimuthwise 
plots  of  Figure  4-7,  the  vector  plots,  with  the  exception  of  a  few 
isolated  points,  do  not  show  this  same  degree  of  marked  change.  Speci¬ 
fically,  the  two  per  rev  signal,  which  is  quite  pronounced  at  point  516 
in  the  X  and  Z  directions  (Figure  4-7)  due  to  the  close  proximity  of 
the  probe  to  the  rotor  tip  path  plane,  undergoes  only  gradual  changes 
when  viewed  in  the  X-Z  plane  of  vectors.  For  instance,  the  change  in 
Vx  from  a  blade  angle  of  ip  =  50°  to  ip  =  100°  is  about  -10  ft/sec.  The 
change  in  for  this  same  interval  is  -5  ft/sec.  In  vector  form  the 
magnitude  of  the  resultant  changes  from  44  ft/sec  to  49  ft/sec  and  the 
angle  change  is  about  15°.  With  the  vector  shown  every  5°  this  change 
takes  place  over  a  span  of  10  pictures  between  =  50°  and  ip  -  100°. 

The  difference  between  the  two  different  representations  is  further 

amplified  by  the  expanded  scale  of  the  V  ,  V  ,  V  plots  of  Figure  4-7 

x  y  2 

as  opposed  to  the  relatively  small  scale  of  the  vector  plots  of  Figures 

4-1  to  4-6.  Therefore  the  azimuth  varying  plots  of  V  ,  V  and  V  pre- 

x  y  z 

sent  the  greatest  amount  of  information  about  the  instantaneous  veloci¬ 
ties. 
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APPENDIX  A 

FOURIER  SERIES  AVERAGING  METHOD 

Early  in  this  work  many  problems  were  encountered  in  attempts  to 
reproduce  the  data  of  previous  tests.  A  possible  solution  originated 
from  the  suspicion  that  the  data  reduction  technique  of  the  computer 
program,  COSPX,  was  in  error.  Specifically,  it  was  believed  that  the 
method  of  arithmetically  averaging  the  voltage  signals  from  each  sensor 
over  14  revolutions  using  the  assumption  of  periodicity  may  have  been 
incorrect.  This  would  be  the  case  if  the  period  of  the  signal  was  not 
a  constant  within  reasonable  limits  (~  +  5°)  because  a  cancellation 
effect  would  result  from  the  averaging  process. 

In  an  effort  to  eliminate  this  susoected  problem,  a  Fourier  Series 
averaging  technique  was  considered.  In  principle  this  approach  entails 
fitting  a  Fourier  Series  to  the  signal  for  each  of  the  14  revolutions 
of  the  rotor  in  order  that  the  waves  may  be  shifted  prior  to  their 
addition  and  subsequent  averaging.  The  analysis  for  this  method 
proceeds  as  follows. 

In  order  to  obtain  the  Fourier  coefficients,  a  system  subroutine 
is  used  which  accepts  as  input  the  data  points  and  returns  the  coeffi¬ 
cients  a„  and  b  of  equation  A-l. 
n  n 

F ( t ) i  =  aQ  +  an  cos(^-  t)  +  bn  sin(^  t)  A-l 

If  the  number  of  data  points  per  revolution  is  N  then  the  number  n  is 
equal  to  N/2  as  stipulated  by  the  subroutine.  So,  in  this  particular 
case  n  =  35. 
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Equation  A-l  can  also  be  expressed  in  the  form 

F(t)i  s  fo  +  Cn  sin(?T  1  +  *n}  A_2 

where  — s - j 

Cn  an  +bn 

0n  =  tan"1  | 

Recall  that  there  are  14  equations  in  the  form  of  A-2,  one  for  each 
revol ution. 

Next,  in  order  to  match  the  peaks,  an  average  value  of  <f>n  for  the 
first  harmonic  is  obtained  from 


14 

-  hi 

♦i '  Mr-  « 

where  *  the  average  first  harmonic  phase  shift 

Once  the  ^  value  is  calculated,  each  wave  can  be  phase  shifted  to 
$1  by  shifting  all  terms  by  the  amount  A<f>.  given  by 


A$.j  =  -  4>i  i  =  1 ,14 

Now  each  wave  can  be  reconstructed  in  the  following  manner 


A-4 


F(t) 


i shi fted 


=  C. 


36 

Z 

n«l 


t  +  4, 


A^) 


where  n  =  1 ,36 

The  final  step  is  to  reduce  the  14  equations  to  1  by  simply  arithmeti 
cally  averaging  the  values  of  CR  and  (<f>n  +  A<j>. ). 


C 


1 


14 

I 

i=l 


14 


A- 5 


For  example: 
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Certainly  all  the  first  terms  will  have  the  sane  value  of  <+>n  since  this 
was  dictated  earlier  by  equations  A-3  and  A-4,  but  the  remaining  35  <>n 
values  are  averaged  as  in  equation  A-5. 

In  principle  the  method  appears  valid  owing  to  the  fact  that  the 
fundamental  frequency  in  the  signal  is  the  prominent  2  per  revolution 
blade  frequency.  Though  in  itself,  the  method  did  not  give  good 
results,  it  spurred  interest  in  obtaining  knowledge  of  the  intermediate 
results  in  the  program  COSPX.FTN.  This  interest  led  to  the  development 
of  a  program  which  plots  the  voltage  signal  for  each  of  the  14  revolu¬ 
tions  against  blade  angle.  Shown  in  Figure  A-l  are  the  voltage  traces 
of  two  revolutions  from  sensor  A  and  the  average  obtained  from  all  14 
revolutions.  This  plot  shows  only  a  small  cancellation  effect  if  any, 
and  consequently  the  method  was  deemed  adequate.  This  is  not  to  say 
that  the  Fourier  Series  method  may  not  provide  even  better  results 
provided  the  believed  software  problem  can  be  eliminated.  One  inherent 
problem  exists  in  using  this  approach  and  that  is  the  limitation  of  the 
system  frequency  response  to  the  number  of  harmonics  times  the  blade 
frequency  which  is  -  1300  Hz.  It  was  in  part  for  this  reason  that 
the  compensators  mentioned  in  2.2  were  removed  from  the  system. 


Figure  A-l.  Arithmetic  Averaging  of  Blade  Revolutions 
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APPENDIX  B 

REDESIGNED  ANEMOMETER  BRIDGE  CIRCUIT 

Redesign  of  the  DISA  55D05/102C  circuits  was  prompted  by  circuit 
transients  which  resulted  in  probe  destruction  on  two  occasions.  An 
operational  circuit  for  use  as  an  alternative  to  these  units  is  pre¬ 
sented  in  Figure  B-l.  The  circuit  is  essentially  the  same  as  the  DISA 
circuits  except  that  the  individual  transistors  have  been  replaced  by 
an  RCA  CA3140AE  operational  amplifier  and  an  ARCHER  VN67AF  power 
transistor.  The  power  transistor  is  necessary  to  supply  the  compen¬ 
sating  current  to  the  bridge. 

Although  this  circuit  is  notably  simplified  its  performance  was 
not  sufficiently  improved  to  warrant  replacement  of  the  DISA  circuits. 
For  instance,  though  frequency  response  was  about  equal,  the  new 
circuit  measured  2  to  5  times  the  noise  level  (light  bulb  tested  in 
place  of  probe)  when  measured  using  the  Spectral  Dynamics  UA500 
spectrum  analyzer.  Due  to  lack  of  time  the  DISA  units  were  repaired 
and  perfection  of  the  new  circuit  was  not  pursued.  One  suggestion  for 
improvement  of  the  circuit  would  be  the  substitution  of  a  Teledyne 
model  1426  op-amp  which  has  a  lower  noise  specification  than  the  RCA 


op-amp. 


Redesigned  Constant  Temperature  Anemometer  Control  Circuit 
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APPENDIX  C 
COMPUTER  PROGRAMS 

The  programs  presented  here  deal  with  the  various  aspects  of  data 
manipulation  necessary  to  produce  the  vector  plots  and  azimuth-varying 
velocity  plots  of  Figs.  4-1  to  4-7.  The  function  of  each  program  and 
important  points  are  discussed  at  the  beginning  of  each  program  list¬ 
ing.  The  listings  are  in  alphabetical  order. 

The  program  EZPLOT.FTN  was  used  to  create  the  plots  of  Fig.  4.7. 
The  tabular  data  for  these  curves,  which  is  also  provided  by  the 
program,  is  given  in  Fig.  C-l. 

The  curves  of  Fig.  A-l  were  created  by  running  VOLTS. FTN  to  plot 
the  average  of  the  14  revolutions  and  EZVOL.FTN  to  plot  each  individual 
revolution. 

The  programs  SHIFT. FTN  and  VECTOR. FTN  were  used  to  plot  Figs.  4-1 
through  4-6. 
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Figure  C-l  continued 
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Figure  C-i  concluded 
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THIS  PROGRAM  IS  DESIGNED  TO  MANIPULATE  THE  DATA  OP  AN 
EXISTING  INSTANTANEOUS  DATA  SET  WITH  THE  END  RESULT  BEING 
AN  AVERAGE  DATA  SET  CALLED  "AVCASE.  DAT**.  THE  PROGRAM 
PIRST  DECODES  THE  DATA  OP  “CASESH.DAT"  WHICH  IS  IN  ITS 
9  DIOIT  INTEGER  FORM  IT  THEN  ADDS  ALL  THE  VX. VY  AND  VZ 
DATA  POR  THE  REVOLUTION  AND  DIVIDES  BV  THE  NUMBER  OP  DATA 
POINTS  IN  THE  REVOLUTION.  (CURRENTLY  SET  TO  EVERY  10  DEGREES) 
THE  DO  40  STATEMENT  CAN  BE  CHANOED  IP  THE  USER  WISHES  SOME 
OTHER  DEGREE  INCREMENT.  POR  INSTANCE:  POR  EVERY  5  DECREES 
THE  STATEMENT  WOULD  READ.  DO  40  AZ-1,72.  1.  ALSO  THE  DIVISORS 
IN  THE  AVERAGING  STATEMENTS  MUST  BE  CHANCED 

ONCE  THE  AVERAGE  VX.  VY. VZ  HAVE  BEEN  DETERMINED  THE  RESULT 
POR  EACH  POINT  IS  CODED  ONCE  AGAIN  IN  THE  9-DIGIT  INTEGER 
FORM  IN  THE  DATA  FILE  "AVCASE.  DAT". 

IN  ORDER  TO  USE  THIS  PROGRAM.  IT  MUST  FIRST  BE  COMPILED 
AND  LINKED.  THE  PROCEDURE  FOR  THIS  IS  AS  FOLLOWS: 

YOU  COM  AVG 
COMPUTER  RETURNS  > 

YOU  FTB 

COMPUTER  RETURNS  FTB> 

YOU  AVQ-AVQ 
COMPUTER  RETURNS  FTB> 

YOU  / 

COMPUTER  RETURNS  ENTER  OPTIONS 
YOU  MAXBUF-4800 
COMPUTER  RETURNS  FTB> 

YOU  // 

COMPUTER  RETURNS  > 

YOU  RUN  AVG 

THE  PROGRAM  WILL  NOW  RUN  AND  THE  DO  LOOP  VARIABLES 
WILL  APPEAR  ON  THE  SCREEN  TO  INDICATE  THE  STATUS. 


INTEGERS  IDATAC 72) •  AZ.  1X1.  IYt.  IZt.  1X2,  IY2,  IZ2.  19 
1  . NDATA< 1188),  IPT 

0PEN(UNIT-2» NAME* 'ARRAY.  DAT'.  TYPE* 'OLD ACCESS- 'DIRECT ' ) 
OPEN (UNIT-3,  NAME* 'AVCASE.  DAT '.  TYPE- 'NEW ' . ACCESS- 'DIRECT ' 
1  . RECORDS! ZE-1 ) 


CALL  INI TT < 960) 
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DO  100  IPT-1.1188 
WRITE <  5,  DIPT 

1  F0RMAT(5X,  'POINT  «  ',I4> 

1X2-0 

IY2-0 

IZ2-0 

DO  40  AZ-2,  72,  2 

READ<2 ' IPT) < IDATA<N) , N=l» 72) 

I 9* I DATA (AZ) 

I X 1-19/ 1000000 

I Y 1 * ( I 9- 1X1*1 000000 ) / 1 000 

I Z1-I9-IX 1*1000000-1 Yl*1000 

IX1-IX1+IX2 

IY1-IY1+IY2 

IZ1-IZ1+IZ2 

IX2-IX1 

IY2-IY1 

IZ2-IZ1 

40  CONTINUE 

C  WR  ITE  ( 5,  * )  I XI ,  IY1,  IZ1 

C  ***  REMEMBER  TO  CHANGE  THE  DIVISOR  IN  THE  FOLLOWING 
C  ***  THREE  STATEMENTS  WHEN  CHANGING  THE  NUMBER  OF 
C  ***  POINTS. 

IX1-IX2/36 

IY1-IY2/36 

IZ1-IZ2/36 

C  WRITE<5. *>IX1, IY1. IZ1 

NDATA< IPT )«1000000*IX 1+1000*1 Yl+IZi 
WRITE(3' IPT)NDATA< IPT) 

WR ITE ( 5. 35)NDATA< IPT) 

35  FORMAT ( 3X,  'NDATA  -  ',110) 

100  CONTINUE 

CL0SE(UNIT-1 > 

CLOSE< UNIT-2) 

CLOSE (UN IT-3) 

STOP 

END 
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COSPX.FTN 

C  ***************************************************** 

THE  PROGRAM  COSPX. FTN  IS  DESIGNED  TO  GATHER  INSTANT 
VELOCITY  DATA  FOR  AN  ENTIRE  CASE.  THE  OPERATOR 
INSTRUCTIONS  ARE  THOSE  USED  WHEN  TAKING  AN  EN¬ 
TIRE  CASE.  THIS  DOES  NOT  MEAN  THAT  THE  PROGRAM 
CANNOT  BE  USED  TO  TAKE  DATA  FOR  OTHER  THAN  AN 
ENTIRE  CASE.  IF  LESS  DATA  IS  BEING  TAKEN  JUST 
IGNORE  THE  INSTRUCTIONS  FOR  LOCATING  THE  PROBE. 

THE  DATA  REDUCTION  METHOD  USED  IS  THE  COSINE 
TO  THE  P  POWER  METHOD.  THE  PROGRAM  CAN  BE  SET 
TO  USE  EITHER  CONSTANT  VALUES  OF  P  FOR  EACH 
SENSOR  OR  TO  USE  AN  ITERATIVE  SCHEME  BETWEEN 
U  AND  P  IF  THE  CALIBRATION  DATA  WARRANTS  THIS 

THE  RESULT  OF  THE  PROGRAM  IS  A  DATA  SET  CALLED 
ARRAY. DAT  WHICH  IS  A  DIRECT  ACCESS  FILE  AND 
THEREFORE  CANNOT  BE  PRINTED  OUT  USING  THE  LIST 
COMMAND.  THE  STORED  DATA  IS  IN  A  9-DIGIT  INTEGER 
CODE.  WHEN  DECODED*  BY  PROGRAMS  SUCH  AS  E2PLOT  OR 
VECTOR  >THE  RESULT  IS  VX,VY  AND  VZ  FOR  EACH  POINT 
IN  THE  TUNNEL.  THERE  ARE  11BB  RECORDS  IN  THE 
DATA  SET  ARRAY.  DAT.  EACH  RECORD  HAS  A  RECORDS  I ZE 
OF  72  SUCH  THAT  VX, VY  AND  VZ  CAN  BE  RECORDED  AT 
5  DEGREE  INCREMENTS  OF  THE  ROTOR  FOR  ONE  ROTOR 
REVOLUTION. 

C  ******************************************************* 
c  *****  THE  FOLLOWING  VARIABLES  ARE  DECLARED  INTEGER*4 
C  *****  BECAUSE  THEY  ARE  INVOLVED  IN  THE  9-DIGIT  INTEGER 
C  *****  CODING  SCHEME. 

INTEGER*4  ARRAY < 72 >, NLOC, NVEL, IVX, IVY; IVZ 
1  ,  I;  IX,  I  Y,  IZ 

C  *****  THE  ARRAY  BUFFER  IS  USED  TO  STORE  THE  VOLTAGES 
C  *****  DELIVERED  BY  THE  A/D  CONVERTER  FROM  CHANNELS  A 
C  *****  B  AND  C  OF  THE  ANEMOMETER. 

INTEGER  BUFFER ( 1024, 3) 

C  *****  THE  ARRAYS  EAANC,  EBANC  AND  ECANG  ARE  USED  TO 

C  *****  STORE  THE  VOLTAGES  FROM  A  B  AND  C  FOR  14  ROTOR 

C  *****  REVOLUTIONS. 

C  *****  THE  ARRAY  V  WILL  BE  USED  TO  STORE  VEFF  A,  VEFF  B 

C  *****  AND  VEFF  C.  THE  ARRAY  U  WILL  BE  USED  TO  STORE  THE 

C  *****  VALUES  OF  UA. UB,  AND  UC. 


DIMENSION  EAANC < 14 ) , EDANG < 14 ) .  ECANG ( 1 4 > ,  V<3).  U<3> 

1  , CSINE<3).  SNE<3) 

WRITE (  5.  11  ) 

II  FORMAT < 3X,  '  IS  THIS  THE  BEGINNING  OF  A  CASE?  ',/, 

1  '  IF  YES  PRINT  A  Y  .  IF  NOT  PRINT  A  N  '  ) 

READ<5, 111)  ANS 

III  FORMAT  < 1 A1 ) 

IF  (ANS.  EQ.  '  Y  ' )  GO  TO  1112 

QPEN<UNIT=1,  NAME*  'OLDLDC.  DAT  ',  TYPE= 'OLD', ACCESS*  'DIRECT  '  ) 
0PEN(UNIT=2,  NAME*  'ARRAY.  DAT  ',  TYPE* 'OLD ' ,  ACCESS*  'DIRECT  '  ) 
OPEN! UNI T=3,  NAME*  'NLQC.  DAT',  TYPE* 'OLD', ACCESS* 'DIRECT  '  > 
READ<  1  '  1  )  I,  IX,  IY,  12,  IYCNTS 
GO  TO  1123 
C 

1112  I X=28 

I  Y=20 
I  Z  =  12 
1*1 

I YCNTS*0 

OPEN(UNIT=l,  NAME* 'OLDLDC.  DAT', TYPE* 'NEW', ACCESS* 'DIRECT  ' 

1  , FORM* 'UNFORMATTED  ' , RECORDSI ZE*5  > 

0PEN(UNIT*2,  NAME* 'ARRAY.  DAT', TYPE* 'NEW ' , ACCESS*  'DIRECT  ' 

1  ,  FORM* 'UNFORMATTED', RECORDSI ZE=72> 

OPEN ( UNI T=3, NAME*  'NLOC.  DAT ' , TYPE* 'NEW',  ACCESS*  'DIRECT  ' 

1  ,  FORM* 'UNFORMATTED', RECORDSI ZE=1 ) 

C 

1123  IF(IY  GE.  12)  WRITE(5,  113)  IX,  IY,  IZ 

113  FORMAT <  3X,  '  THE  PROBE  SHOULD  BE  NOW  ON  HOLDER  A',/, 

1  '  AT  LOCATION:  ',  //,  '  X*  ',  13,  /.  '  Z=  ',  13,  /,  '  Y=  ',  13) 

C 

IF<  IY.  LE.  8.  AND.  IY.  GE.  4)  WRITE<5,  114)  IX,  IY,  IZ 

114  FORMAT ( 3X.  '  THE  PROBE  SHOULD  BE  NOW  ON  HOLDER  B',/, 

1  'AT  LOCATION:  ',  //,  '  X*  ',  13.  /.  '  Y*  ',  13,  /,  '  Z=  ',  13) 

C 

IF(  IY.  LE.  0.  AND.  IY.  GE  -8)  WRITE(5,  115)  IX,  IY,  IZ 

115  FORMAT < 3X,  '  THE  PROBE  SHOULD  BE  NOW  ON  HOLDER  C',/. 

1  '  AT  LOCATION:  ',  //,  '  X*  ',  13,  /,  '  Y=  ',  13,  /,  '  Z=  ',  13) 

C 

IF<  IY.  LE.  -12)  WRITE(5,  116)  IX,  IY,  IZ 

116  FORMAT ( 3X, '  THE  PROBE  SHOULD  BE  NOW  GN  HOLDER  D',/, 

1  '  AT  LOCATION:  ',  //,  '  X*  ',  13,  /,  '  Y=  '.  13,  /,  '  Z*  ',  13) 
WR ITE <  5, 650) 

650  FORMAT < 3X,  '  WHEN  IT  IS  THERE  PRINT  OK') 

READ<5, 651)  LOC 

IF (LOC. NE.  'OK')  GO  TO  1123 

651  FORMAT  < 1A2) 


O-  u  c 


MX»(2B-IX)/4+l 
MY* (20-1 Y ) /4+1 
MZ*( 12-1 Z ) /3+1 

C  *****  THE  SUBROUTINE  SET ADI  IS  USED  TO  SET  THE  PARAMETERS 
C  *****  OF  THE  A/D  CONVERTER. 

CALL  SET  ADI  (3.  3,  1.  5) 

DO  101  JY=MY, 11 
DO  100  JZ=MZ,  6 
DO  99  JX=MX, 18 
C 

.1  ND*1 

IF  ( ABS<  I X  > .  LE.  12.  AND.  IZ.  LE.  3.  AND.  ABS(IY)  LE  12) 

C  *****  THE  SUBROUTINE  ZERO  SETS  ALL  POINTS  IN  THE  ROTOR 
C  *****  INTERFERENCE  VOLUME  EQUAL  TO  ZERO. 

1  CALL  ZERG<ARRAY. IND) 

IF<  IND.  EQ.  0)  GO  TO  91 
652  WRITE<5, 117) 

117  FORMAT <3X, '  NOW  PUSH  THE  SWITCH  ON  INTERFACE  ') 

C  *****  THE  SUBROUTINE  GETAD  IS  USED  TO  ACTUATE  THE  A/D 
C  *****  CONVERTER  TO  FILL  THE  ARRAY,  BUFFER 


CALL  GETAD (BUFFER) 

WRITE <5, 5000) 

5000  FORMAT < 3X, '  BUFFER  WAS  FILLED') 

DO  85  J=l, 72 
C  (72*360/5) 

DO  9  M=l.  3 
I  A= J 

DO  8  JU=1 , 14 

C  14-REVOLUTIONS 

IF(M.  EQ.  1)  EAANG< JO >=FLOAT< BUFFER < IA,  1 >>/ <32768  /B  > 
IF(M.  EQ.  2)  EBANG(JJ) “FLOAT ( BUFFER ( IA,  2) ) /( 32768.  /8  ) 
IF (M  EQ.  3)  ECANG  <  JU ) “FLOAT ( BUFFER ( I A,  3) > / (32768.  /S.  ) 
IA*IA+72 
8  CONTINUE 


C  *****  THE  SUBROUTINE  AVERAG  IS  USED  TO  GET  AVERAGE 
C  *****  VOLTAGES  FOR  ONE  REVOLUTION  FROM  THE  14  REVS 
C  *****  STORED  PREVIOUSLY. 


IF (M.  EQ.  1) 
IF(M.  EQ.  2) 
IF(M.  EQ.  3) 
CONTINUE 
VOLTAGE  TO 


CALL  AVERAG ( EAANG,  AVA) 
CALL  AVERAG < EBANG. AVB > 
CALL  AVER AC ( ECANG. AVC ) 

VELOCITY  CALCULATION 
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USE  OF  CALIBRATION  CURVES  TO  OBTAIN  THE  CORRESPONDING 
EFFECTIVE  VELOCITIES 

CALIBRATION  CONSTANTS  FOR  PARTICULAR  PROBE 


A0=  7552 
Al-i.  7513 
A2-.  7304 
A3=.  01033 
A4=  001053 

V<  1  )*A0+A1*AVA+A2*AVA**2+A3*AVA**3*A4*AVA**4 

BO*.  1367 
B  1  *  1 . 833 
B2=.  8454 
B3*~*.  02726 
B4*. 004034 

V ( 2 ) *BO+B 1 *AVB+B2*AVB**2+B3*AVB**3+B4*AVB**4 

CO*.  6043 
Cl*l.  8375 
C2=.  7809 
C3=-  01365 
C  4*.  003839 

V ( 3 ) =CO-*-C  1*AVC+C2*AVC**2+C3*AVC**3+C4*AVC  **4 

WRITE  STATEMENTS  TO  PRINT  VOLTAGES  AND  UA  UB  UC 
*****  NOTICE  THAT  THESE  STATEMENTS  HAVE  A  "D"  IN  COLUMN 
*****  i.  THIS  MEANS  THAT  IF  WHEN  COMPILING  THE  USER  TYPES 
*****  ••  COM  COSPX  /DE  "  (DON'T  INCLUDE  THE  QUOTES)  THE  "D" 
*****  STATEMENTS  WILL  BE  INCLUDED  IN  THE  PROGRAM.  IF  /DE 
*****  IS  NOT  USED  THE  PROGRAM  TREATS  D  STATEMENTS  AS 
*****  COMMENT  STATEMENTS. 

IF ( J.  GT.  7 1 )  WR ITE  < 6»  288>AVA,  AVB,  AVC 
288  FORMAT (2X<  'AVA  =  ',F7.  3,  '  AVB  *  '.F7.  3,  '  AVC  *  ',F7.  3) 
IF(  J.  QT.  71 ) WRITE  (6.  289 )V(  1),  V<2>.  V(3) 

289  FORMAT <3X,  'VI  =  ',F7.  3,  '  V2*  ',F7.  3,  '  V3=  ',F7.  3) 

1F(  J.  GT.  71 )  WRITE  (  5,  290JAVA.  AVB,  AVC 

290  F0RMAT<2X,  'VOLTAGES  -  3F7.  3) 

1 F  ( J.  GT.  71  )WRITE  <  5,  291  )  V<  1  ) ,  V < 2 ) ,  V( 3 ) 

291  FORMAT <4X.  'VI  =  ',F7.  3,  '  V2«  ',F7.  3,  '  V3*  ',F7.  3) 


CALCULATE  MAGNITUDE  OF  THE  ABSOLUTE  VELOCITY 
*****  NOTICE  THAT  THE  CALL  STATEMENT  IS  CURRENTLY  COMMENTED 
*****  OUT  OF  THE  PROGRAM  BUT  IF  AN  ITERATION  BETWEEN  U  AND 
*****  p  WERE  NECESSARY  THE  STATEMENT  CAN  BE  EASILY  CHANGED. 


oooo  *-nonoooo  a  ooo 
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CALL  UP ITER (V.  VEL<  P»  J) 

P*.  78 

VEL  IS  THE  STREAM  VECTOR  -U- 
VEL*8  <  (V(  1  )**<2.  /P  )  +V  <  2 )  **  <  2.  /P)*V<3)**(2.  /P 
1  )  )  /2.  >**<P/2.  > 

IF<  J.  GT.  71 )  WRITE  <6,  287>VEL 
IF(J.  GT  71 ) WRITE (5,  287)VEL 
FORMAT <  3X.  'VEL  =  '.F7.  3) 

FIND  ANCLES  IN  PROBE  SYSTEM  OF  COORDINATES 
THE  VARIABLE  K  REPRESENTS  A  B  AND  C. 

DO  15  K= 1 ,  3 
I F ( K.  EQ.  1  >P=.  76 
IF ( K.  EQ.  2)P=.  80 
I F ( K.  EQ  3>P*.  78 
IF(  J.  GT.  71 ) WRITE <6,  296 )P 
296  FORMAT ( 3X,  '  P  =  ',  F5.  3) 

CSINE<K)»<V(K)/VEL>**<1.  /P) 

IF(CSINE<K) .  GT.  1.0)  CSINE<K)  =  1.0 
IF<CSINE(K).  GT.  0.  99.  AND.  J.  GT.  71  ) WRITE (5,  297) 

297  FORMAT (3X<  'COSINE  IS  GREATER  THAN  .99') 
ACSN«AC0S(CSINE(K>  > 

SNE  <  K ) =S I N  <  ACSN ) 

IF <  J.  CT.  71  )WRITE<5,  295)SNE(K) 

I F  ( J.  CT.  7 1 )  WRITE  <  6.  295)SNE(K) 

295  FORMAT <3X,  'SINE  =  ' ,  F6.  3) 

IF <  J. CT.  71 ) WRITE (6,  286>CSINE(K> 

296  F0RMAT<5X.  'CSINE(K)  =  ',F6.  3) 

FIND  PERPENDICULAR  COMPONENTS  OF  VELOCITY 

IN  PROBE  SYSTEM  OF  COORDINATES 
5  U<K.)«VEL*SNE<K> 

IF<  J.  CT.  71 )  WRITE  (6,  292 ) U <  1  ) ,  U< 2 ) ,  U<  3 ) 

292  FORMAT <3X<  'U1  »  ',F7.  3,  '  U2  *  ',F7.  3,  '  U3  *  ',F7.  3) 
IF(  J.  CT.  71 ) WRITE  <  5.  293 >U<  1 ),  U<2),  U(3> 

293  FORMAT <3X*  'U1  *  ',F7.  3,  'U2»  '.F7.  3.  'U3-  '-F7.3) 

IN  TUNNEL  SYSTEM  OF  COORDINATES 
VZ=0. 57735*  <  U ( 1 ) +U ( 2 ) -U ( 3 ) ) 

VY*  7071*(U< 1 )-U<2> ) 

VX»-0.  40825*(U( 1 >+U<2) )-0.  8165*U(3> 


C 

C  ***** 

D 
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C  CORRECTION  FACTORS  INCORPORATED  IN  VY  AND  VZ 
C  AS  FUNCTIONS  OF  VEL. 

I F( ABS( VY ) .  LT.  1.  0)G0  TO  299 
VYCOR=VEL*(-.  091  )+.  58 
VY=VY+VYCOR 

299  IF(ABS(VZ>.  LT.  1.  0)G0  TO  300 
VZCOR=VEL*  < .  108)-.  06 
VZ=VZ+VZCOR 

C  *****  THE  FOLLOWING  STATEMENTS  ARE  THE  9-DIGIT  INTEGER 
C  *****  CODING  SCHEME  USED  TO  CONSERVE  MEMORY  SPACE. 

300  I VX=IFI X  <  VX*5.  +0.  5>+499 
I  VY= IFI X  ( VY*5.  +0.  5 )  +499 
IVZ=IFIX  <  VZ*5.  +0.  5 )  +499 
NVEL= I VX * 1000000+ IVY* 1000+ I VZ 
ARRAY ( J ) =NVEL 

D  IF<  J.  GT.  60)WRITE<6,  69>VX,  VY,  VZ 

69  FORMAT < 3X«  'VX='»  IX,  F7.  2,  2X,  'VY*',  IX,  F7  2,  2X, 

1  'VZ=',  IX,  F7.  2) 

IFI  J.  CT.  71 )  WRITE  (5,  68)VX,  VY,  VZ 
68  FORMAT  <  3X,  'VX  ='.  F7.  2,  2X,  'VY  =',  F7.  2,  2X,  'VZ  =  ',  F7  2) 

85  CONTINUE 

91  NLOC=< I X +499 ) * 1 000000+ ( IY+499)*1000+I Z+499 
C 

WR ITE ( 2 ' I )  (ARRAY (K) , K=l. 72) 

WR I TE  <  3  '  I )  NLOC 

IF (  IND.  EQ.  0)  GO  TO  990 

WRITE  <  5,  991) 

991  FORMAT ( 3X ,  '  WAS  THE  PROBE  IN  THE  CORRECT  LOCATION'' ',  //, 

1  'IF  YES  PUSH  RETURN  .  IF  NOT,  ENTER  -N-  ') 

READ <5, 111)  BACK 
IF<BACK.  EQ.  'N '  >  CO  TO  652 
990  1=1+1 

I X= I X-4 

I F (  IX.  EQ.  12.  AND.  IZ.  LE.  3.  AND.  ABS<  IY>.  LE.  12)  IND=0 
IFdND.  EQ.O)  WRITE<5,  86) 

86  FORMAT ( 3X,  '  YOU  ARE  NOW  IN  THE  DANGEROUS  AREA',/, 

1  '  MOVE  THE  PROBE  TO  Z=6  (78  CNTS) , AND  THEN  TO  X=-16'> 
IF(  IX.  LT. -40)  IX-28 
IFdND.  EQ.  1)  WRITE!  5,  96  >  IX 
96  FORMAT (3X,  '  MOVE  PROBE  TO  LOCATION',//, 

1  '  X=',  14) 

I F (  IX.  EQ.  28.  AND.  IZ.  LE  3  AND.  ABSIIY).  LE.  12)  WRITE<5,  106) 
106  FORMAT <3X,  '  BEFORE  YOU  MOVE  THE  PROBE  TO  X*28',/, 

1  '  TAKE  IT  HORIZONTALY  TO  Z«6  (78  CNTS)') 

WRITE! 5, 92) 

92  FORMAT (3X»  '  IF  YOU  WISH  TO  CONTINUE  PUSH  RETURN',/, 

1  '  IF  NOT  PRINT  AN  — N — .  ') 

READ<  5, 111 )  ANS 
IF(ANS.  EQ.  'N')  GO  TO  1000 

99  CONTINUE 
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MX-1 

IZ-IZ-3 

IF<  IZ.  LT.  -3)  IZ=12 
IZCNTS«< 12-1 Z) *13 
WRITE<5, 124)  IZ, IZCNTS 

124  FORMAT ( 3X, '  MOVE  PROBE  HORIZONTALY  TO',//, 

1  '  Z-  '.14.'  EQUIVALENT  TO',  14,  '  (CNTS)') 

125  WRITE<5,  126) 

126  FORMAT ( 3X , '  PUSH  RETURN  TO  CONTINUE') 

READ(5,  111)  ZANS 

IF  ( ZANS.  EQ.  'N'>  00  TO  125 

100  CONTINUE 
MZ  =  1 

I Y= I Y-4 

I YCNTS=I YCNTS+52 

IF<  IY.  EQ.  8.  OR.  IY.  EQ.  0.  OR.  IY.  EQ.  -12)  IYCNTS*0 
350  FORMAT <3X>  '  MOVE  THE  PROBE  VERTICALLY  TO  ',  /, 

1  '  Y*  ',14.'  EQUIVALENT  TO',  14,  '  (CNTS)') 

IF <  IY.  EQ.  8 )  WRITE<5,  360) 

360  FORMATOX,  '  CHANGE  PROBE  TO  HOLDER  — - - - 

IF(IY.  EQ.  0)  WRITE<5.  370) 

370  FORMATOX.  '  CHANGE  PROBE  TO  HOLDER  — - - - 

IF(  IY.  EQ.  -12)  WR ITE (  5,  300) 

380  FORMATOX,  '  CHANGE  PROBE  TO  HOLDER  — - - - 

WRITE(5,  350)  IY,  IYCNTS 
390  WRITE <5,  400) 

400  FORMATOX.  '  PRINT  RETURN  TO  CONTINUE  ',/, 

l  '  YOU  CAN  NOT  TERMINATE  THE  SESSION  FROM  HERE') 

READ <5, 111)  YANS 

IFCYANS.  EQ.  'N')  CO  TO  390 

101  CONTINUE 

1000  I ZCHEK=I Z-3 
I XCHEK=I X— 4 

IF ( I XCHEK.  LT.  -40.  AND.  IZCHEK.  LT.  -3)IZ«12 
IF( IXCHEK.  LT.  -40.  AND.  IZCHEK.  LT.  -3)  IY=IY-4 
WRITE(l'l)  I, IX, IY, IZ, IYCNTS 
DO  1001  IU=1,3 

1001  C  LOSE  <  UN I T= I U ) 

C  LOSE ( UN I T*6  > 

STOP 

END 

C  SUBROUTINE  ITERATES  BETWEEN  U  &  P  TO  GET  VEL=U 

C 

C  SUBROUTINE  UPITER<V, VEL, P, J) 


oo  r>  nonooooooDonnnnonooooorjoorj 
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DIMENSION  V(3)< VELI (16) 

FIRST  GUESS  FOR  P  IS  0.  6 
P».  6 

VEL*< <V< 1 )**<2.  /P)+V<2)**<2.  /P)+V<3)**<2  /P) )/2. 

1  )**<P/2.  ) 

VELI < 1 )=VEL 
DO  77  1-1.15 
RE  =  <  VEL-.  346)/.  885 
P*. 593+  0042*RE 

VELI < 1  +  1 )  =  < <V< 1 )**<2.  /P  >  +V  ( 2 )  **  <  2.  /P  >  +V ( 3 ) ** 

1  <2.  /P>  )/2.  )  **  ( P/2.  ) 

IF  <  ABS<  VELI  <  1  +  1  )  -VELI  <  I  )  ).  GT.  0.  01.  AND.  I  EQ  15  AND. 
1  <U.  GT.  71  )  >WRITE<5.  283) 

283  FORMAT ( 3X<  '***  P  HAS  NOT  CONVERGED') 

IF < ABS< VELI < 1  +  1 ) -VELI <I)).  LT.  0.  01)  GO  TO  78 

77  CONTINUE 

78  VEL-VELI ( I > 

I F ( J.  GT  71 )WRITE<5. 410>P 
410  FORMATOX,  'P  =',F7  4) 

RETURN 

END 

*****  NOTE  THE  FOLLOWING  FUNCTIONS  ARE  USED  IN  THE  K**2 
*****  METHOD  AND  HAVE  BEEN  COMMENTED  OUT 

FUNCTIONS  FOR  K=F<U>  FOR  EACH  SENSOR 

SENSOR  A 
FUNCTION  UKKU) 

AUO=  2796 
AU1 .  03397 
AU2*2.  387E-3 
AU3*-6.  541E-5 
AU4»5.  959E-7 

UK 1 »AU4*U**4+AU3*U**3+AU2*U**2+AU 1 *U+AUO 

RETURN 

END 

SENSOR  0 

FUNCTION  UK25U) 

BUO*.  06048 
BUI— 7.  794E-4 
BU2-5.  249E-4 
BUS— 2.  231E-5 
BU4-2.  441E-7 

UK2«BU4*U**4+BU3*U**3+BU2*U**2+BU1*U+BU0 

RETURN 

END 


oooo  or>  r>  o 
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SENSOR  C 

FUNCTION  UK3(U) 

CUO— .  2898 
CU1*.  06627 
CU2*-3.  443E-3 
CU3-7.  156E-5 
CU4=-5. 227E-7 

UK3=CU4*U**4+CU3*U**3+CU2*U**2«-CU1*U+CU0 

RETURN 

END 


SUBROUTINE  AVERAG< A,  AV) 

DIMENSION  A< 14),  SUM< 14) 

SUM ( l )*A< 1 ) 

DO  9  1*1,  13 

9  SUM( I  +  l )*SUM< I )+A< 1  +  1  ) 

AV=SUM<  1 4 )  / 14. 

RETURN 

END 

SUBROUTINE  TO  SET  VX*VY=VZ=0 

WHEN  THE  PROBE  IS  IN  THE  DANGEROUS  REGION 

SUBROUTINE  ZERO < ARRAY, IND) 

INTEGER*4  ARRAY (72), NVEL,  1VX,  IVY, 3VZ 
DO  10  1*1.  72 
I VX=499 
1 VY=499 
IVZ*499 

NVEL*IVX*1000000+IVY*1000+IVZ 
10  ARRAY ( I )*NVEL 

IND*0 

,  RETURN 
END 


on  on  on  on 
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EZPLOT.FTN 


C  THE  PROGRAM  EZPLOT.FTN  IS  DESIGNED  TO  PLOT  THE  AZIMUTH 
C  VARYING  VELOCITIES  FOR  X,  Y  AND  Z.  THE  RESULT  OF  THE 
C  PROGRAM  IS  NOT  A  PLOT  BUT  A  FILE  IN  THE  DIRECTORY  CALLED 
C  F0R001.DAT  WHICH  WILL  CONTAIN  VX,  VY  OR  VZ  DEPENDING  ON 
C  WHICH  YOU  SPECIFY  IN  THIS  PROGRAM  WHEN  YOU  RUN  IT 

C  THE  PROGRAM  WILL  ASK  YOU  FOR  A  PARTICULAR  POINT  NUMBER 
C  WHICH  MUST  BE  FROM  1  TO  1188.  IF  YOU  NEED  HELP  IN  FINDING 
C  THE  POINT  NUMBER  OF  A  TUNNEL  COORDINATE*  SEE  REFTAB. FTN 

C  IN  ORDER  TO  GRAPH  THE  FILE  FOROOl.DAT  YOU  MUST  FOLLOW  THIS 
C  THIS  SIMPLE  PROCEDURE  TO  USE  EASY  GRAPHING. 

1)  AFTER  THE  PROGRAM  HAS  BEEN  RUN  TYPE  THE  FOLLOWING 
"COPY  FOROOl.DAT  PT.  DAT"  (OMIT  THE  QUOTES) 

2)  AFTER  THE  PROMPT  HAS  BEEN  RETURNED  TYPE  IN 
"EZC" 

3)  THE  SCREEN  WILL  BE  CLEARED  AND  A  BELL  SOUNDED  THEN  TYPE 
"RUN  VEL" 

4)  IF  YOU  WANT  TO  CHANGE  THE  SCALES  OR  ADD  TITLES  SEE 
THE  MANUAL  ON  EASY  GRAPHING 

C  *********************************************************** 

INTEGER*4  ARRAY (72)*  CODVEL*  IVX.  IVY,  IVZ 
DIMENSION  ANGL(72*  3) 

0PEN(UNIT*2.  NAME- 'ARRAY.  DAT '*  TYPE-  'OLD ' , ACCESS-  'DIRECT ' ) 

1  CALL  INITT (960) 

WRITE(  5.  2) 

2  FORMAT (3X,  '  PLEASE  ENTER  "I"  A  DATA  POINT') 

READ ( 5*  * )  I 

C  ***  THE  FOLLOWING  IS  A  DIRECT  ACCESS  READ  STATEMENT 
C  ***  TO  READ  THE  DATA  FROM  ARRAY  OF  UNIT  2  CORRES- 
C  ***  PONDING  TO  THE  POINT  NUMBER  THE  USER  HAS  SPECIFIED. 

READ (2  'I )  (ARRAY< J), J-l* 72) 
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DO  10  11*1,72 

IF  ( 1 1.  LT.  2) GO  TO  201 

GO  TO  202 

201  WRITE<5,  99) 

99  FORMATOX,  'DO  YOU  WISH  TO  INCLUDE  A  TUNNEL  CORRECTION 
1  FACTOR  ?  //3X,  'TYPE  “Y"  FOR  YES  OR  "N"  FOR  NO') 

READ (5, 199) REPLY 
199  FORMAT (1A1) 

IFCREPLY.  EG.  'N')GO  TO  202 

WRITE ( 5» 103) 

103  FORMATOX,  'TYPE  IN  Y-CORRECTION  ' ) 

READ< 5, 110) YCOR 

110  FORMAT  <F3.  1 ) 

WRITE ( 5,  104) 

104  FORMATOX,  'TYPE  IN  Z-CORRECTION '  ) 

READ (5, 120) ZCOR 

120  FORMAT  <F3.  1 ) 

202  CODVEL*ARRAY  < 1 1 ) 

C  ***  THE  FOLLOWING  STATEMENTS  ARE  THE  DECODING  ROUTINE 
C  ***  FOR  THE  9-DICIT  INTEGERS. 

VX*  <  CODVEL/ 1000000-499 > /5 
I  VX*  <  CODVEL/ 1 000000 ) * 1 000000 
VY*( FLOAT ( (CODVEL-IVX )/1000)-499.  >/5 
I  VY*( ( CODVEL-IVX ) / 1000) *1000 
V  Z -FLOAT  <  CODVEL- 1 VX - 1 VY-499 ) / 5. 

8  ANGL< II,  1 )*VX 

VY-VY+YCOR 
VZ-VZ+ZCOR 
ANGL( 1 1,  2)=VY 
ANGL  <11,  3)*VZ 
10  CONTINUE 

WRITE<5,  55)  ( ANGL  <  M,  1),  M»l.  72,  1 ) 

55  FORMATOX,  6F9.  2) 

WRITE15, 56) 

56  FORMAT <  1 X,  / ) 

WRITE  (  5,  57)  <  ANGL  <M,  2),  M-l,  72,  1 ) 

57  FORMATOX,  6F9.  2) 

WRITE (5, 58) 

58  FORMAT < IX.  / ) 

WRITE (  5,  59 )  <  ANGL  < M,  3 ) ,  M*  1 ,  72,  1) 

59  FORMATOX,  6F9.  2) 

WRITE< 5,  400) 

400  FORMATOX,  'IF  YOU  WANT  TO  PLOT  VX,  VY  OR  VZ  PUSH  RETURN' 

1  ,  /,  '  IF  NOT  TYPE  "N"  ') 

READ (5, 401) KEY 

401  FORMAT ( 1 A1 ) 
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IF  (KEY.  EQ.  'N')GOTO  1000 
WRITE<5. 500) 

500  FORMAT (3X,  'TYPE  IN  VX, VY.  OR  VZ  '  > 

READ(  5,  501 )  ANS 

501  FORMAT <1A2) 

I F ( ANS.  EQ.  'VX  ' )  J-l 
IF  (ANS.  EQ.  'VY')  J-2 
IF  (ANS.  EQ.  'VZ')  J*3 
IF(ANS.  EQ.  'VX')  GOTO  777 
IF  (ANS.  EQ.  'VY')  GOTO  888 
DO  100  1*1 ,  72 
YZ*ANGL< I.  J) 

WRITE(1,  *>YZ 

100  CONTINUE 

GOTO  999 

777  DO  101  11*1/  72 

IX«( ( 1 1-1 >*5>*2 
YX*ANGL( II.  J) 

WRITE ( 1. *)YX 

101  CONTINUE 

GOTO  999 

888  DO  102  12*1,72 

YY*ANGL(I2.  J) 

WRITE( 1.  *>YY 

102  CONTINUE 

999  CONTINUE 

1000  WRITE(5.  901) 

901  FORMAT (3X. 'WOULD  YOU  LIKE  TO  PLOT  A  DIFFERENT  POINT? 
1  TYPE  "Y"  FOR  YES  AND  "N"  FOR  NO'  ) 

READ (5, 902) LET 

902  FORMAT ( 1 A1 ) 

IF  (LET.  EQ.  'Y'>  GOTO  1 
CLOSE (UN IT-2) 

STOP 


EZVOL.FTN 


C  ***************************************** 

C  THE  PROGRAM  EZVOL.  FTN  CREATES  A  FILE  IN 
C  THE  DIRECTORY  CALLED  F0R003.  DAT  WHICH 
C  CONTAINS  THE  VOLTAGES  FOR  EACH  REV  OF 
C  EITHER  CHANNEL  A  B  OR  C.  PLOTTING  CAN 
C  BE  DONE  USING  THE  EASY  GRAPHING  ROUTINE. 

C  THE  PROGRAM  WHEN  RUN  WILL  ASK  IF  YOU  WANT 
C  THE  VOLTAGES  FROM  A,  B  OR  C,  AND  OF  WHICH 
C  REVOLUTION.  IF  YOU  WANT  TO  PLOT  THE  AVER- 
C  AGE  OF  ALL  THE  REVS  YOU  MUST  RUN  THE  PRO- 
C  GRAM  VOLTS.  FTN. 

C  BEFORE  YOU  CAN  PLOT  THESE  VOLTAGES  YOU  MUST 
C  RUN  THE  PROGRAM  MAKEBUF. FTN  TO  CREATE  THE 
C  DIRECT  ACCESS  FILE  BUFFER.  DAT. 

C  ****************************************** 

I NTECER  BUFFER  < 1 024, 3 > 

DIMENSION  FNT<72, 14,3) 

0PEN<UNIT-1.  NAME-  'BUFFER.  DAT',  TYPE® 'OLD ' ,  ACCESS= ' 

1  DIRECT' > 

READ(l'l)  <BUFFER<N<  1)»N=1,  1024) 

READ< 1 '2)  <BUFFER<N< 2), N-l, 1024) 

READ<1'3)  <BUFFER<N< 3). N=l,  1024) 

DO  5  M=1 ,  3 
K-0 

DO  10  J-l.  14 
DO  20  1*1,72 
L-K+I 

C  ***  DECODING  THE  BINARY  VOLTAGE  VALUE. 

FNT <  I ,  J,  M ) -FLOAT  <  BUFFER  <  L,  M)  )  /  < 32768.  /8.  ) 

20  CONTINUE 

K-K+72 

10  CONTINUE 

5  CONTINUE 

C  ***  INITT <960)  CLEARS  THE  PAGE  AND  SETS  BAUD 

CALL  INITT < 960) 

WRITE <5,  500) 

500  FORMAT <4X.  'DO  YOU  WANT  TO  PLOT  VOLTAGES  A,  B  OR  C  .  ' 
READ<  5, 501 ) REPLY 

501  FORMAT  < 1A1 ) 

IF <REPLY.  EG.  'A')  GOTO  100 
IF < REPLY.  EQ.  'B'>  GOTO  200 
IF < REPLY.  EQ.  'CM  GOTO  300 
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C  THIS  IS  THE  STORAGE  ROUTINE  FOR  CHANNEL  -C- 

300  WRITE <5, 502) 

502  FORMAT <4X. 'WHICH  REVOLUTION  WOULD  YOU  LIKE  TO  PLOT. 

1  1,2,3.  .  .  14  ?') 

READ  <5,  * )  I  REV 
JC*IREV 

DO  301  10=1, 72 
YC  =  ( FNT ( I C , JC , 3 ) ) 

WRITE (3,  * ) YC 

301  CONTINUE 
GOTO  700 

C  THIS  IS  THE  STORAGE  ROUTINE  FOR  CHANNEL  -B- 

200  WRITE <5, 511) 

511  FORMAT ( 4X, 'WHICH  REVOLUTION  WOULD  YOU  LIKE  TO  PLOT. 

1  1,  2,  .  .  .  14?') 

READ <5, * ) IREVB 
JB*IREVB 

DO  201  IB=1 »  72 
YB=FNT  < IB. JB, 2) 

WRITE <3,  *>YB 

201  CONTINUE 
GOTO  700 

C  THIS  IS  THE  STORAGE  ROUTINE  FOR  CHANNEL  -A- 

100  WR1TE(5, 512) 

512  FORMAT ( 4X, 'WHICH  REVOLUTION  WOULD  YOU  LIKE  TO  PLOT 
1  ,  1,2...  14?') 

READ (5, *) IREVA 
JA=IREVA 

DO  101  I A=1 , 72 
YA*FNT  ( I  A,  JA,  1) 

WR ITE ( 3,  *) YA 

101  CONTINUE 

700  CLOSE (UNIT=1 ) 

STOP 

END 
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LOCK . FTU 

C  ******* *************** ******* ******** ************ 

C  THE  PROGRAM  LOOK.  FTN  ALLOWS  THE  USER  TO  INSPECT  THE 
C  9-DIGIT  INTEGER  DATA  STORED  IN  THE  PILE  CASESH. DAT 
C  THIS  WOULD  PRIMARILY  BE  USED  IN  A  DIAGNOSTIC  SENSE. 

C  THE  FILE  CASESH.  DAT  IS  THE  RESULT  OF  THE  RUNNING  0 
C  SHIFT  FTN.  SINCE  THIS  REPRESENTS  AN  ENORMOUS  AMOUNT 
C  OF  DATA  ONLY  ONE  AZIMUTH  POSITION  IS  ASKED  FOR  AT 
C  ONE  TIME.  THE  DATA  FOR  THIS  AZIMUTH  WILL  BE  STORED 
C  IN  THE  FILE  F0R006  DAT. 

C  **************************************************** 

INTEGER*4  I DATA < IS. 6. 1 1 ), ARRAY ( 1 14 ), AZ 

OPEN<UNIT*l .  NAME® 'CASESH  DAT TYPE® 'OLD ' ,  ACCESS® 'DIRECT ' ) 
WRITE(  5.  100) 

100  FORMATOX,  'TYPE  IN  THE  AZIMUTH  OF  THE  DATA(I£:  1  TO  72>'> 
READ  <5.  *)  AZ 

READ< 1 'AZ) ( <  < IDATA(K.  I »  J  > ,  K* 1 ,  18).  I  * 1 .  6 ) ,  J® 1 .  11) 

WRITE (6. 200) I DAT A 
200  FORMATOX.  9110) 

CLOSE < UNI T=1  ) 

STOP 

END 
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MAKE3UF .  FT*J 

C  *******##***##•###*****##**•»■*#****##*#****** 

C  THE  PROGRAM  MAKEBUF.  FTN  IS  DESIGNED  TO  FILL 
C  THE  DIRECT  ACCESS  FILE  BUFFER.  DAT  WITH  THE 
C  VOLTAGES  FROM  A,  B.  AND  C  FOR  ALL  14  REVS. 

C  THE  PROGRAM  MUST  BE  USED  WITH  A  TRIGGER  TO 
C  ACTIVATE  THE  A/D  CONVERTER.  THE  WIND  TUNNEL 
C  INTERFACE  UNIT  OR  A  WAVETEK  MAY  BE  USED. 

C  TO  SUBSEGUENTLY  PLOT  THE  DATA  GATHERED  HERE. 

C  RUN  EITHER  VOOLTS.  FTN  FOR  THE  AVERAGE  OR 
C  E2VOL.  FTN  FOR  EACH  INDIVIDUAL  REV. 

C  REMEMBER  TO  USE  THE  FTB  LINK  PROCEDURE  AND 
C  USE  A  MAXBUF=4752. 

C  *******************•»*«***'»*»**«***'********** 

INTEGER  BUFFER < 1024, 3) 

0PEN(UNIT=1,  NAME® 'BUFFER.  DAT', TYPE=  'NEW ' , ACCESS®  'DIRECT  ' 
1  , FORM® 'UNFORMATTED  '» RECORDSI ZE=1024 ) 

C  ***  SET ADI  READIES  THE  A/D  CONVERTER 

CALL  SETAD1  (3.  3,  1.5) 

WRITE<  5,  200) 

200  FORMATOX,  'NOW  PUSH  THE  SWITCH  ON  THE  INTERFACE') 

C  ***  GETAD  ACTUALLY  FILLS  BUFFER 

CALL  GETAD (BUFFER) 

WR ITE ( 5, 300) 

300  FORMATOX,  'BUFFER  HAS  BEEN  FILLED  ) 

C  ***  BUFFER  IS  BEING  WRITTEN  INTO  UNIT  ONE 
WRITE< I  '1 )  (BUFFER  <N,  1 ), N=l,  1024) 

WR ITE  < 1  '2 )  < BUFFER <N, 2). N-l.  1024) 

WRITE< 1  '3 >  (BUFFER(N, 3),  N=l,  1024) 

READ ( 1  '  1 )  < BUFFER ( N,  1 ). N=l,  1024) 

READ ( 1  '2 )  <BUFFER(N. 2), N-l.  1024) 

READd'3)  < BUFFER  <N,  3).  N-l.  1024) 

C  ***  THE  WRITE  STATEMENT  PRINTS  EVERY  100TH  VALUE 
C  ***  OF  BUFFER  ON  THE  SCREEN. 

WRITE  (  5.  400)  <  <  BUFFER  <  M,  MN),  M-l,  1024,  100),  MN*1, 3) 

400  F0RMAT(/1X, 5110, /,  51 10) 

CLOSE (UN IT® l ) 

STOP 

END 
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REFTAB. FTN 

C  #*•**#**##***##*#*#*#*•**#**#*■**##*****#•*********** 

C  THE  PROGRAM  REFTAB.  FTN  CREATES  A  TABLE  IN  THE 
C  DIRECTORY  UNDER  THE  NAME  F0R002.  DAT.  THIS  TABLE 
C  CONTAINS  THE  X#  Y  AND  Z  COORDINATES ( IN  INCHES)  OF 
C  EACH  OF  THE  1188  POINTS  IN  A  CASE.  THIS  PROGRAM 
C  CAN  BE  LINKED  BY  USING  "LINK  REFTAB"  SINCE  IT  IS  NOT 
C  USING  THE  DIRECT  ACCESS  MODE.  THE  FTN  FILE  MUST 
C  ALSO  BE  COMPILED  USING  "COM  REFTAB"  BEFORE  IT  CAN 
C  BE  RUN.  TO  RUN  THE  PROGRAM  JUST  TYPE  "RUN  REFTAB" 

C  IF  YOU  HAVE  ALREADY  COMPILED  AND  LINKED  THE  PROGRAM 

C  *#***#*###**#*****#**##****#*#*•****#*##**##********# 

DIMENSION  LIST ( 1 188.  4 ) 

1*0 

DO  80  J»l.  11 
DO  70  K*1  #  6 
DO  60  L*1 < 18 
I-I  +  l 

DO  50  M*»l *  4 
WRITEC  5.  *)I,  J»  Kj  L.  M 
IF(M.  EQ.  1  )LIST  <  I  *  M)*I 
IF<M.  EQ.  2 )  L I  ST  ( I  <  M )  =28-  <  L- 1 ) *4 
IF<M.  EQ.  3>LIST<I,M>=20-(J-1>*4 
IF(M.  EQ.  4>LIST(I,M)»12-(K-1)*3 
50  CONTINUE 
60  CONTINUE 
70  CONTINUE 
80  CONTINUE 

WRITEC2.  90) < (LIST ( I.  M) .  M=1 , 4 ) ,  1  =  1 ,  1 188 ) 

90  FORMATdOX,  4110) 

STOP 

END 
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REWRIT. FTN 


THIS  PROGRAM  ALLOWS  DATA  TO  BE  WRITTEN  OVER 
POINTS  THAT  HAVE  ALREADY  BEEN  TAKEN.  THE 
USER  MUST  SPECIFY  THE  POINT.  I,  THE  POSITION 
IX,  IY,  12  AND  THE  VERTICAL  COUNTER  READING, 
IYCNTS.  THEN  INSTVEL  CAN  BE  RUN  AS  USUAL 
USING  THE  “OLD  CASE"  ANSWER. 


***  THE  VALUES  IX, IY. IZ  ARE  IN  INCHES 
INTECER*4  I,  IX.  IY,  IZ 

OPEN<UNIT«l,  NAME* 'QLDLOC.  DAT',  TYPE® 'NEW ',  ACCESS- 'DIRECT ' 
1  . FORM- 'UNFORMATTED ', RECORDSIZE-5) 

WRITE<5,  10) 

10  FORMATOX,  'FIRST  TYPE  IN  "I",  THE  PT.  #  YOU  WISH  TO 
1  REENTER') 

READ<5.  *)I 
WRITE( 5, 20) 

20  FORMAT <3X,  'NOW  TYOE  IN  THE  CORRESPONDING  X, Y.  AND  Z 
1  ',/,  'POSITIONS  IN  INCHES.  ') 

READ  <5.  *)  ANS 

IF(ANS.  EG.  'EMGOTO  100 

READ  (  5,  *)IX.  IY,  IZ 

30  WRITE<5.  40) 

40  FORMATOX,  'NOW  TYPE  IN  THE  -Y-  COUNTER  READING 
1  FOR  YOUR  POSITION') 

READ< 5,  *) IYCNTS 

WRITE(1 '1)1,  IX,  IY.  IZ,  IYCNTS 

CL0SE<UNIT*1 ) 

100  STOP 
END 
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SHIFT. FTN 


.  .  .  THE  PROGRAM  SHIFT  REARRANGES  THE  STORED  FORMAT  OF  THE 
.  .  .  DATA  IN  THE  FILE  ARRAY.  DAT  CREATED  BY  RUNNING  THE 
.  ..  PROGRAM  COSPX.  FTN.  INSTVEL.  FTN.  VOLTS.  FTN  OR  AVERAGE. 

.  .  .  FTN.  FOR  SUBSEQUENT  USE  IN  THE  PROGRAM  VECTOR.  FTN 

.  .  THE  COUNTER  THAT  APPEARS  ON  THE  SCREEN  AS  SHIFT  IS 
.  .  RUN  GOES  FROM  1  TO  72  AND  CAN  TAKE  CONSIDERABLE  TIME 
.  .  .  TO  ACCOMPLISH  THIS  DEPENDING  UPON  THE  NUMBER  OF  POINTS 
.  .  BEING  SHIFTED.  (IE:  AN  X-PLANE  HAS  1 1  X  6  OR  66  POINTS 
.  .  .  WHEREAS  AN  ENTIRE  CASE  HAS  11  X  6  X  IS  OR  1168  POINTS) 

•a**#*##***###*#*##*#**#*#**####******#*****###***#*###***** 


INTEGER*4  IDATA(18, 6.  1 1 ). ARRAY ( 1 14 ) .  I TEST ( 18.  6.  1 1 > 

BYTE  ID1 (80).  ID2(80) 

0  WRITE(5.  100) 

00  FORMAT ( '  TYPE  IN  SIX  LETTER  SOURCE  DAT  FILE  NAME: ',*> 
READ(  5.  101.  ERR*  10.  END-1000)  I,  ( ID1  (  J).  J-4,  1+4) 

101  FORMAT ( Q. 12A1) 

.  .  .  UNLESS  THE  UNIT  NUMBER  IN  THE  DATA  TAKING  PROGRAM 
.  .  .  HAS  BEEN  CHANGED.  THE  UNIT  NUMBER  REQUESTED  HERE 
.  .  WILL  BE  "2". 

WRITE ( 5.  102) 

102  FORMAT ( '  UNIT  NUMBER  WHICH  SOURCE  FILE  IS  IN:  *> 

READ(5. Ill)  IU1 

WRITE (5. 106) 

READ(5. 107)  NIVEL 
121  WRITE<5. 104) 

104  FORMAT ( '  MIN/MAX  9-POINTS  ALONG  X-Y-Z  AXES  (IE:  1.18, 

1  1.6.1.11)  '*) 

***  THE  "IE"  ABOVE  IS  FOR  AN  ENTIRE  CASE 

READ(  5.  *.  ERR-121.  END-1000)  IXS.  IXE.  IYS,  I  YE.  IZS.  IZE 

.  .  .  THE  ARRAY  LENGTH  ASKED  FOR  HERE  WILL  BE  "72"  UNLESS 

.  .  .  THE  DATA  TAKING  PROGRAM  IS  RECORDING  MORE  THAN  72 

.  .  .  POINTS  PER  BLADE  REVOLUTION  OR  UNLESS  YOU  ARE  SHIFTING 

.  .  .  AVERAGE  DATA  IN  WHICH  CASE  THE  ARRAY  SIZE  IS  1. 

106  FORMAT ( *  ARRAY  LENGTH  OF  SINGLE  DATA  SET:  *> 

107  FORMAT (115) 

11  WRITE<5, 108) 

108  FORMAT ('TYPE  IN  SIX  LETTER  DESTINATION  .DAT  FILE  NAME: 
1  *) 

READ ( 5.  109.  ERR-1 1. END-1000)  I,  ( ID2( J) . J-4,  1+4) 

FORMAT (Q,  12A1 ) 
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.  .  .  THE  UNIT  NUMBER  ASKED  FOR  HERE  CAN  BE  ANY  UNIT  OTHER  THAN 
.  .  .  THAT  BEING  USED  BY  THE  SOURCE  FILE  AND  NOT  5,6  OR  7. 
WRITE<5,  110) 

110  FORMAT < '  UNIT  NUMER  WHICH  DESTINATION  FILE  IS  TO  BE 
1  STORED  ',♦) 

READ! 5. Ill)  IU2 

111  FORMAT  < 1 1 1 ) 

ID1 < 1 )*'S ' 

ID2( 1 )*  'S ' 

ID1 (2)* 'Y ' 

I D2 ( 2 ) *  '  Y  ' 

ID1 (3)* ' : ' 

I D2 ( 3 )  *  ' :  ' 

IDKI+4)*'.  ' 

I D2 ( I +4 ) * ' .  ' 

ID1  < 1+5)" 'D ' 

ID2( I+5>* 'D' 

ID1 ( I+6>* ' A ' 

I D2 ( I +6 )  *  '  A ' 

ID1 ( 1+7)*  'T  ' 

ID2(I+7)*'T' 

NPTS-0 

...  THESE  THREE  DO  LOOPS  JUST  PERFORM  THE  MULTIPLICATION  OF 
...  THE  MAXIMUM  X, Y, 2  LIMITS  TO  GIVE  THE  NUMBER  OF  POINTS 
.  .  .  BEING  SHIFTED. 

DO  12  I-IZS.  I2E 
DO  13  K-IYS.  I YE 
DO  14  J-IXS.  IXE 
NPTS-NPTS+1 
14  CONTINUE 

13  CONTINUE 

12  CONTINUE 

WRITE<5. 199)  NPTS 
199  FORMAT ( '  #-POINTS*'» 110. ♦) 

OPEN < UNIT* I U2* NAME-ID2. TYPE* 'NEW '» ACCESS* 'DIRECT ' 

1  .  FORM* 'UNFORMATTED ' »  RECORDSI ZE-NPTS ) 

OPEN ( UN I T- 1 U 1 « NAME-ID1.  TYPE- 'OLD'.  ACCESS* 'DIRECT ' ) 

DO  1  11*1. NIVEL 
111*0 

DO  2  K-12S.  I2E 
DO  3  I-IYS,  I YE 
DO  4  J-IXS,  IXE 
III-III+l 

READ<IU1 'III)  <ARRAY<JJ). JJ*1.  NIVEL) 

I  DATA  (J.  I.  K) -ARRAY <  II) 

4  CONTINUE 

3  CONTINUE 

2  CONTINUE 
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300 

1 


301 

C 

1000 


WRITE< IU2 'll)  < ( ( I DATA (K, I.  J>,  K*IXS,  IXE) .  I=IYS,  I YE), 
1  J-IZS,  IZE) 

WRITE  (  5,  300)  ARRAY  <  II ),  IDATAdXE,  I  YE,  IZE),  II 
FORMAT<  IX,  110,  2X,  110.  3X.  'I  *  ',12) 

CONTINUE 

THE  FOLLOWING  WRITE  STATEMENT  IS  COMMENTED  HERE 
BUT  CAN  BE  USED  BY  REMOVING  THE  "C"  TO  PRINT  OUT 
THE  ENTIRE  IDATA  ARRAY  WHICH  CONTAINS  THE  9-DIGIT 
CODED  FORM  OF  THE  X, Y, Z  VELOCITIES. 

WRITE <6,  301) IDATA 
FORMAT <3X,  6110) 

CONTINUE 

CLOSE (UN IT* I U1 > 

CL0SE(UNIT*IU2) 

END 


STORE. FTN 

C  *******#*#*#*#*****#*##*#***############•»#*# 

C  THE  PROCRAM  STORE.  FTN  IS  A  SIMPLE  READ  WRITE 
C  OPERATION  WHICH  SPLITS  A  WHOLE  CASE  OF  DATA 
C  INTO  TWO  SEPARATE  FILES  IN  ORDER  TO  ALLOW 
C  STORAGE  ON  TWO  FLOPPY  DISKS.  A  CASE  OF  DATA 
C  OSES  669  BLOCKS  OF  STORAGE  AND  A  FLOPPY  DISK 
C  HOLDS  ONLY  ABOUT  430. 

C  THE  RESULT  OF  THE  PROGRAM  IS  TWO  DATA  FILES 
C  CALLED  HALF1.DAT  AND  HALF2.DAT  WHICH  CAN  BE 
C  STORED  SEPARATELY  ON  TWO  FLOPPIES.  TO  RE- 
C  CREATE  THE  FILE  ARRAY.  DAT  FOR  USE  IN  PLOTTING 
C  THE  FOLLOWING  PROCEDURE  SHOULD  BE  FOLLOWED. 

C  TYPE  IN  "MERGE  HALF2.DAT  HALF1.DAT” 

C  THIS  CREATES  A  FILE  HALFi.DAT  WHICH  IS  THE 
C  ADDITION  OF  HALF2  TO  HALF 1 . 

C  NEXT  JUST  TYPE  COPY  HALF!.  DAT  ARRAY.  DAT. 

C  ***###**#####*####*#*##*#*********##********* 

INTEGER *4  ARRAY<72),I 

0PEN(UNIT«2,  NAME* 'ARRAY.  DAT'.  TYPE* 'OLD ' ,  ACCESS 
1  DIRECT') 

0PEN(UNIT*3»  NAME* 'HALF2.  DAT',  TYPE* 'NEW '•  ACCESS 
1  DIRECT',  FORM* 'UNFORMATTED',  RECORDS I ZE*72> 
OPEN! UNI T*6,  NAME* 'HALF1 .  DAT ',  TYPE* 'NEW ' ,  ACCESS 
1  DIRECT ' .  FORM* ' UNFORMATTED ' .  RECORDSI ZE*72 ) 

DO  100  1*1,763 

READ (2 ' I )  <  ARRAY  <  J ) , J* 1 , 72 ) 

WRITE (6' 1 )  < ARRAY( J) ,  J-1,72) 

100  CONTINUE 

DO  101  1*764,  USB 

R£AD<2 ' I )  <ARRAY< J), J*l, 72) 

11*1-763 

WRITE(3 ' 1 I )  (ARRAY! J) , J*l, 72) 

101  CONTINUE 
CLOSE (UN IT-2) 

CL0SE(UNIT*3) 

CLOSE (UN IT-6) 

STOP 

END 
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VECTOR. FTN 

C  ************************************************** 

C  THE  PROGRAM  VECTOR.  FTN  ALLOWS  THE  USER  TO  PLOT 
C  THE  INSTANTANEOUS  VELOCITY  VECTORS  IN  THE  X.  Y 
C  OR  Z  PLANES.  IN  ORDER  TO  USE  THIS  PROGRAM  THE 
C  DATA  AS  TAKEN  IN  THE  TUNNEL  MUST  BE  SHIFTED  USING 
C  THE  PROGRAM  SHIFT. FTN.  THIS  REARRANGES  THE  DATA 
C  IN  A  DIRECT  ACCESS  FILE  WHICH  HAS  ALL  THE  VELOCITY 
C  DATA  FOR  A  SPECIFIC  BLADE  ANGLE  INSTEAD  OF  ALL  THE 
C  ANGLES  FOR  ONE  POINT. 

C  *************************************************** 

C  TO  RUN  VECTOR.  FTN  IT  IS  FIRST  NECESSARY  THAT  THE 
C  PROGRAM  BE  COMPILED  AND  LINKED  SINCE  THE  DATA 
C  IN  THE  PROGRAM  IS  DIRECT  ACCESS  THE  ORDINARY  “LINK1' 
C  COMMAND  CANNOT  BE  USED.  THE  FTB  OR  TKB  DIRECT 
C  ACCESS  LINK  PROCEDURES  MUST  BE  USED.  THIS  IS  DONE 
C  AS  FOLLOWS: 

C  YOU  FTB 
C  COMPUTER  FTB> 

C  YOU  VECTOR-VECTOR 
C  COMPUTER  FTB> 

C  YOU  / 

C  COMPUTER  ENTER  OPTIONS: 

C  COMPUTER  FTB> 

C  YOU  MAXBUF-4800 
C  COMPUTER  FTB> 

C  YOU  // 

C  *****  NOTE  THAT  THE  BUFFER  HAS  BEEN  DECLARED  AS  4800 
C  *****  THIS  IS  THE  SIZE  NEEDED  WHEN  ACCESSING  AN 
C  *****  ENTIRE  CASE  OF  DATA.  THE  BUFFER  SIZE  MUST  BE 
C  *****  AT  LEAST  4  TIMES  THE  RECORDSIZE  OF  THE  FILE 
C  *****  4  TIMES  1188  EQUALS  4752 

C  **************************************************** 

INTEGER*4  IDATA1 (18, 11). IDATA2(18, 11) 

LOGICAL  ITEST, IHRD 

COMMON  /Cl/IHRD,  ITEST, IU, INVELS, INVELE,  INVEL 
COMMON  /C6/IXP, IYP, IHOLD 
I HOLD- 1 
ITEST-.  TRUE. 

C  ***  THE  SUBROUTINES  ATTACH,  INITT,  TERM,  CHRSIZ 
C  ***  ,  ANMARM,  GRLOOK,  NUMBUF,  NUMDRW,  AND  KEYFUN 
C  ***  ARE  ALL  PLOT  10  SYSTEM  SUBROUTINES 
CALL  ATTACH 
CALL  INITT < 960 ) 

CALL  TERM <2, 1024) 

CALL  CHRSIZ <4 ) 
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C 

CALL  ANMARM 
CALL  6RL00K ( 1 ) 

CALL  NUMBUF (2) 

CALL  NUMDRW<3) 

CALL  KEYFUN< 1 ) 

C 

CALL  ANMODE 
CALL  ASK1< IFILE) 

IF  (  IFILE.  EQ.  2 )  WR  I TE  <  5.  15) 

SHIFT.  FTN  IS  A  PROGRAM  WHICH  REARRANGES  THE  DATA 
TO  AN  X,Y ,1  AT  A  PARTICULAR  AZIMUTH  FOR  EACH  POINT 
RATHER  THAN  THE  AS  TAKEN  FORM  OF  X,  Y,  Z  FOR  ALL 
AZIMUTH  POSITIONS  AT  EACH  POINT. 

15  FORMAT < '  TYPE  IN  "RUN  SHIFT"') 

IF<  IFILE.  EQ.  2)  GO  TO  1000 

20  CONTINUE 

CALL  FILE 
CALL  ASK2 
CALL  NEWPAG 
CALL  BORDR 
CALL  TUNNEL 
I HQLD=2 

CALL  CZAXIS<0) 

CALL  HOME 
CALL  ANMODE 
CALL  WHERE 
CALL  SETPAR 
CALL  INFO IN 
CALL  STPANM 
IBUF-1 
CALL  NEWPAG 
C 

UJ**INVELS 
GO  TO  91 

19  IF(IHRD)  JJ=JJ+INVEL 

IF(  (JJ.  GT.  INVELE).  OR.  <  JU.  EQ.  INVELE)  )  JJ=1NVELS 
91  IF(IHRD)  CALL  NEWPAG 

CALL  CZAXIS(O) 

CALL  INFOUT 
CALL  CZAXlS(l) 

CALL  ORIGIN(JJ) 

IF(IHRD)  GO  TO  2002 
C  ***  INVELS-STARTING  BLADE  ANGLE 
C  ***  1NVELE-ENDING  BLADE  ANGLE 
C  ***  INVEL-BLADE  ANGLE  INCREMENT 
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DO  2000  JJ-INVELS. INVELE. INVEL 
C 

2002  CONTINUE 

CALL  DATA(  JJ»  I,  J,  IDATA1 *  IDATA2) 

IF(  .  NOT.  IHRD) CALL  CZAXIS(2> 

IF(IHRD)  CALL  CZAXIS(O) 

IF(.  NOT.  IHRD)  CALL  FRAME ( IBUF,  1  ) 

CALL  PL0T2(I, J,  IDATA1,  IDATA2) 

IF (IHRD)  GO  TO  2000 
CALL  ANIMAT (IBUF) 

IBUF«IBUF+1 

IF(  IBUF.  EQ.  3)  IBUF=1 

2000  CONTINUE 
C 

IF (IHRD)  CALL  STPANM 
READ ( 5.  1928)  IQAZ 
1928  FORMAT (1A1) 

CALL  STPANM 
C 

CALL  HOME 

CALL  ANMODE 

IF(IHRD)  WRITE(5, 8991 ) 

8991  FORMAT ('  NEXT  PLOT?  Y-N  ',*) 

READ< 5,  2001 )  IUYT 

IF(  IUYT.  EQ.  '  Y ' )  GO  TO  19 

2001  FORMAT ( 1 A1 ) 

WRITE (5,  9000) 

READ( 5,  9001 )  IDUM 

9000  FORMAT ('  ANOTHER  RUN?:  Y-N  ',*) 

9001  FORMAT ( 1A1 > 

IF(IDUM.  EQ.  'N ' )  GO  TO  1001 
WRITE <5.  21 ) 

21  FORMAT ('  FROM  THE  SAME  SOURCE  FILE?:  Y-N 

READ  (5,  22)  ITES 
22  FORMAT  < 1 A1 ) 

IF ( ITES.  EQ.  '  Y ' )  ITEST=.  FALSE. 

IF(  ITES.  EQ.  'N ' )  ITEST=.  TRUE. 

IF(  ITES.  EQ.  'N' )  CLOSE (UNIT» IU) 

GO  TO  20 

1001  C  LOSE ( UN I T» I U ) 

CALL  STPANM 
1000  CONTINUE 

END 

C  **«*•»•*»»*«****•**#•*«*#********«**** 


c\»  n 
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C  THE  SUBROUTINE  FILE  ASKS  THE  USER  THE 
C  FILE  NAME  OF  THE  DATA  SET  CREATED  BY 
C  THE  PROGRAM  SHIFT.  FTN. 

SUBROUTINE  FILE 
BYTE  I D < 20 ) 

LOGICAL  I  TEST 

COMMON  /Cl/IHRD, ITEST, IU, INUELS,  INVELE. INVEL 
IF( .  NOT.  ITEST)  GO  TO  1001 

17  WRITE*  5,  19) 

19  FORMAT*'  NAME  OF  SIX  LETTER  DAT  DIRECT  ACCESS  SOURCE 
1  FILE.  '.*) 

READ*  5,  18, ERR=17, END=1001 )  I,  ( ID* J),  J=4,  1+4) 

18  FORMAT  *Q,  12A1 ) 

WRITE*  5,  16) 

16  FORMAT*'  UNIT  NUMBER  SOURCE  DAT  FILE  IS  STORED  IN:  ',  *> 
READ* 5,  *)  IU 
ID* 1 )*'S' 

ID<2)*'Y' 

ID*  3) *  ' :  ' 

ID*  I+4)  = '.  ' 

ID* I+5)*'D' 

ID* I +6)— 'A ' 

ID* I+7)='T' 

OPEN ( UN I T» I U,  NAME* ID, TYPE®  ' OLD ' , ACCESS* 'DIRECT') 

C 

1001  RETURN 
END 

C  *********************** ***************** 

SUBROUTINE  INFOIN 
BYTE  ID1 ( 200) 

INTEGER  IDPL* 10) 

COMMON  /C2/ID1, IDPL 
WRITE*  5,  1 ) 

WRITE* 5,  2) 

WRITE (5,  3) 

1  FORMAT*'  TYPE  IN  ALL  INFORMATION  ABOUT  PLOT') 

FORMAT* '  TO  BE  DISPLAYED  WITH  PLOT  ') 

FORMAT*'  LINE  BY  LINE,  AS  MANY  LINES  AS  NEEDED.  ') 

C 

III«1 

ISS‘1 

WRITE* 5,  5) 

5  FORMAT* 'TYPE  IN  INFO.  WHEN  APPEARS') 

WRITE *5,  6) 

6  FORMAT* 'TYPE  IN:  **  TO  TERMINATE  INPUT.  ') 

900  WRITE* 5,  7) 

III-III+l 

7  FORMAT*'  t) 

READ* 3,  8, END-1001 )  I,  <ID1(II),  II-ISS,  (ISS+I-l)) 

8  FORMAT ( Q,  35A1) 
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IF!  ID1  <  ISS).  EG.  .AND.  IDi  <  ISS+1  ) .  EQ  '*')  GO  TO  901 

ISS-ISS+1 
IDPL < 1 1 1 ) -I 
GO  TO  900 

901  IDPL  < 1 )=I I 1-1 
1001  RETURN 
END 

C  *************************************** 

c  ***  IGD-THE  DISTANCE  BETWEEN  GRID  POINTS 
C  ***  I  SC  X— THE  SCREEN  X  MAGNITUDE  FACTOR  FOR  VECTORS 
C  ***  I SC Y-THE  SCREEN  Y  MAGNITUDE  FACTOR  FOR  VECTORS 
C  ***  ISKX-THE  SKEW  FACTOR  IN  SCREEN  X-DIRECT 
C  ***  ISKY-THE  SKEW  FACTOR  IN  SCREEN  Y-DIRECT 
C  ***  I X=X  COORD.  FROM  SCREEN  LOWER  LEFT 
C  ***  I Y=Y  COORD.  FROM  SCREEN  LOWER  LEFT 

C  ***************************** *************** 

SUBROUTINE  SETPAR 

COMMON  /C3/IXA,  I YB.  IL.  IV,  IP 

COMMON  /C5/IGD, ISCX, ISCY, ISKX, ISKY, IX, IY 

COMMON  /C6/IXP,  IYP,  I  HOLD 

COMMON  /Cl/IHRD, ITEST, IU, INVELS, INVELE, INVEL 
WR ITE<  5,  15) 

WRITE  (5,  16) 

20  WRITE(5,  12) 

WRITE!  5,  13) 

WRITE!  5,  14) 

IX=IX-100 
WR  I  TE  (  5,  8 ) 

WRITE!  5,  7) 

WRITE(5,  9) 

WRITE  <5,  15) 

12  FORMAT ( '  IGRID=SI2E  OF  GRID  (TRY  55)') 

15  FORMAT ( '  ') 

16  FORMAT ( '**********  SET  PARAMETERS  ***********') 

9  FORMAT < '  IX,  IY=  RASTER  REFERENCE  POINTS:  ( IE:  100.  100)  '  ) 

13  FORMAT <  '  SC ALE X= SCALING  FACTOR  ON  X-MAGNITUDE  '  ) 

14  FORMAT < '  SC ALEY=SC AL I NG  FACTOR  ON  Y-MAGNITUDE  '  ) 

8  FORMAT < '  SKEWX® AMOUNT  OF  SKEW  SUBTRACTION  ON  "X"  ') 

7  FORMAT < '  SKEWY=AMOUNT  OF  SKEW  SUBTRACTION  ON  "Y"  ') 

WR  I  TE  <  5,  10) 

WR  ITE  (  5,  11)  IGD,  ISCX,  ISCY,  ISKX,  ISKY,  IX,  IY 
11  FORMAT! IX, 717) 

10  FORMAT!'  INPUT:  IGRID,  SCALEX,  SCALEY,  SKEWX, SKEWY,  IX,  IY  ' ) 
READ!  5,  #,  END-1000,  ERR-20)  IGD,  ISCX,  ISCY,  ISKX,  ISKY,  IX,  IY 
IX-IX+100 

1000  CONTINUE 

WRITE<5,  15) 
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C  ***  NOTE:  IF  PLOTTING  OF  AVERAGE  DATA  IS  REQUIRED,  THE  VALUES 
C  ***  IN  THE  FOLLOWING  STATEMENT  WILL  BE  1,1.1  WHEN  THEY  ARE 
C  ***  ASKED  FOR  IN  THE  INSTRUCTIONS. 

WR ITE ( 5, 30)  INVELS,  INVELE,  INVEL 
30  FORMAT ( '  LENGTH  OF  SINGLE  DATA  POINT  ARRAY! 1 . 72, 1 > : ', 

1  313,  IX,  *> 

READ! 5, *, END=1001, ERR  =  1000)  INVELS,  INVELE.  INVEL 
1001  RETURN 
END 

C  **************************************** 

SUBROUTINE  INFOUT 
BYTE  ID1 <200) 

INTEGER  IDPL< 10) 

COMMON  /C2/ID1# IDPL 
CALL  CZAXIS<0) 

ISS=1 

DO  910  K=2,  IDPL < 1  ) 

CALL  MOVABS < 770,  <&20-K*10>> 

CALL  ANMODE 

WRITE  <5,  850)  <  ID1  <  I  ),  I  =  ISS,  <  IDPL  <  K  )  +  ISS- 1  )  ) 

850  FORMAT  <  Q, 35A1 ) 

ISS*ISS+IDPL(K) 

910  CONTINUE 

RETURN 
END 

C  **************************************** 

SUBROUTINE  M0VE3 < I X,  I Y,  I Z ) 

REAL  OUTX,  OUTY.  OUTZ,  TEMP,  RMAT (4,4) 

INTEGER  IX,  IY,  IZ,  IOUTX,  IOUTY 
COMMON  /C6/IXP,  IYP,  IH 
COMMON  /C8/IXHP, IYHP, I 
COMMON  /RMAT/RMAT 
IF<  IH.  EQ.  2)  CO  TO  1234 
IF  <  IH.  EQ.  1)1  XP=  I  XHP 
IF( IH. EQ. 1)IYP-IYHP 
1234  CONTINUE 

X*FLOAT< IX) 

Y=FLOAT< IY) 

Z “FLOAT  < I Z  > 

0UTX«X*RMAT<1,  1)+Y#RMAT<2,  1)+Z#RMAT<3.  1 >+RMAT<4,  1 ) 

C 

OUTY* X*R MAT < 1,  2)+Y*RMAT<2, 2 ) +Z*RMAT < 3,  2 ) +RMAT (4,2) 
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c 

OUTZ«X#RMAT< 1,  3)+Y*RMAT<2,  3>+Z*RMAT<3,  3 > +RMAT < 4, 3 > 
C 

TEMP*OUTX*RMAT (1,4) +OUTY*RMAT < 2< 4 > +OUTZ*RMAT ( 3, 4)  + 
1  RMAT  <4»  4) 

OUT X=OUTX /TEMP 
OUTY*OUTY/TEMP 
I OUTX* IF I X  < OUTX ) +1 XP 
I OUT Y* IF I X  <  OUT Y )  + 1  YP 
CALL  MOVABS< IOUTX, IOUTY) 

RETURN 

END 

C  **************************************** 

SUBROUTINE  DRAW3 (IX- IV, I Z ) 

REAL  OUTX,  OUTY,  OUTZ,  TEMP,  RMAT (4,4) 

INTEGER  IX. IY, IZ, IOUTX, IOUTY 
COMMON  /C6/IXP, I YP, IH 
COMMON  /C8/IXHP,  IYHP,  I 
COMMON  /RMAT/RMAT 
IF<  IH  EQ.  2)  GO  TO  1234 
I F <  IH.  EG.  1)1  XP=  1  XHP 
I F  ( IH.  EG.  1  )  IYP— 1  YHP 
1234  CONTINUE 

X=FLOAT( IX) 

Y=FLOAT ( IY) 

Z*FLOAT< IZ) 

OUTX  =  X*RMAT< 1,  1 ) +Y*RMAT < 2,  1 )+Z*RMAT(3,  1 )+RMAT(4,  1 ) 

C 

OUTY=X*RMAT< 1,  2 ) +Y*RMAT ( 2. 2 ) +Z*RMAT ( 3,  2>+RMAT<4,  2) 
C 

OUTZ=X*RMAT< 1,  3)+Y*RMAT<2,  3  > +Z*RMAT  (  3,  3)+RMAT(4,  3) 
C 

TEMP=OUTX*RMAT  (1,4)  +OUTY*RMAT  <  2,  4  )  +OUTZ*RMAT  <  3,  4  >  + 
1  RMAT (4,  4) 

OUT X=OUTX /TEMP 
OUT Y*OUTY /TEMP 
I OUTX* IF I X ( OUTX )  +  1 XP 
I OUTY*IFI X ( OUTY ) +1 YP 
CALL  DRWABS( IOUTX, IOUTY) 

RETURN 

END 
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c  **************************************** 

SUBROUTINE  INIMAT(PX,  PY,  PZ,  PH E.  RX,  RY,  RZ) 

REAL  ALPHA,  BETA,  GAMMA,  RMAT  (  4,  4),  PS,  PC,  Nl,  N2,  N3 
INTEGER  PX.  PY,  PZ,  PHE,  RX,  RY.  RZ 
INTEGER *4  I TEMP 
COMMON  /C6/IXP,  IYP,  IH 
COMMON  /C8/IXHP,  IYHP.  I 
COMMON  /RMAT/RMAT 
IF(  IH.  EG.  2)  GO  TO  1234 
IF ( IH.  EQ.  1 )  I XP=  I XHP 
IF <  IH.  EG.  1 )  I YP=  I YHP 
1234  CONTINUE 

ITEMP«PX**2+PY**2+PZ**2 
TEMP»FLOAT < I TEMP ) 

TEMP=SGRT(TEMP) 

N 1 “FLOAT ( P  X ) / TEMP 
N2*»FLOAT  ( PY )  /TEMP 
N3*FL0AT ( PZ ) /TEMP 
PS“S IN (FLOAT (PHE ) /180.  0*3.  14159) 

PC*COS( FLOAT < PHE) /ISO.  0*3.  14159) 

RMAT ( 1 . 1 )*Nl*Nl+( 1. -Ni*Nl >*PC 
RMAT ( 1,  2)«N1*N2*( 1.  -PC)+N3*PS 
RMAT (1,3) =N1 *N3* ( 1 . -PC ) -N2*PS 
RMAT ( 1 ,  4)*0.  0 

IF  (RX.  NE.  0)  RMAT  ( 1,  4)»1.  /FLOAT  ( RX  ) 

RMAT (2,  1)«N1*N2*(1.  -PC)-N3*PS 
R MAT (2,2) »N2*N2+ ( 1 .  -N2*N2 )*PC 
RMAT (2.3) =N2*N3* ( 1 .  -PC)+N1*PS 
RMAT  (2,  4)*=0.  0 

IF ( RY.  NE.  0)  RMAT ( 2,  4)  =  1.  /FLOAT ( RY ) 

RMAT (3,  1 )*N1*N3*< 1.  -PC)+N2*PS 
RMAT (3,  2)»N2*N3*( 1.  -PC)-N1*PS 
RMAT ( 3,  3 ) *N3*N3+ ( 1 .  -N3*N3 ) *PC 
RMAT  (3.  4)=0.  0 

IF(RZ.  NE.  0)  RMAT ( 3, 4)«1  0/FL0AT(RZ> 

RMAT <  4,  1  )»0.  0 
RMAT  (4,  2)«0.  0 
RMAT  (4,  3)=0.  0 
RMAT(4,  4)»1.  0 
RETURN 
END 
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C  **************************************** 

c 

SUBROUTINE  BOR DR 
CALL  HOME 
CALL  DRUABS ( 0,  0 ) 

CALL  DRUABS < 1023. 0) 

CALL  DRWABS< 1023, 780) 

CALL  DRUABS (0,  780) 

RETURN 

END 

C  **************************************** 

SUBROUTINE  ASK1  <K> 

COMMON  /Cl/IHRD,  ITEST.  IU,  INVELS,  INVELE,  INVEL 
URITE<  5,  30) 

30  FORMAT < '  HAS  THE  DAT  BEEN  SHIFTED?  "V"  OR  "N" '  *) 
READ<5,  31)  SH 

31  FORMAT  < 1 A1 > 

IF< SH.  EQ  'Y '  )K=1 
I F ( SH.  EQ.  'N '  )K=2 
RETURN 
END 

C  **************************************** 

SUBROUTINE  ASK 2 

COMMON  /Cl/IHRD, ITEST, IU, INVELS. INVElE, INVEL 
UR  I TE  <  5,  33) 

33  FORMAT  (  '  HARDCOPIES  OR  ANIMATION?:  "H"  OR  ‘'A"',*) 
READ(5, 34)  IUAY 

34  FORMAT  < 1 A 1 ) 

IF ( IUAY.  EQ.  'HO  IHRD=.  TRUE. 

IF  ( IUAY.  EQ.  'A ' )  IHRD*.  FALSE. 

RETURN 

END 

C  **************************************** 

C  THE  INSTRUCTIONS  FOR  THIS  PART  ARE  FAIRLY 
C  CLEAR  EXCEPT  FOR  ONE  POINT.  IF  YOU  ARE 
C  ACCESSING  A  UHOLE  CASE  SPECIFY  THE  BEGINNING 
C  AND  ENDING  LOCATIONS  OF  THE  UHOLE  CASE  EVEN 
C  IF  YOU  ARE  ONLY  PLOTTING  ONE  PLANE. 


C  ***  IXS-X  STARTING  LOCATION 
C  ***  IXE-X  ENDING  LOCATION 
C  ***  IYS-Y  STARTING  LOCATION 
C  ***  IYE-Y  ENDING  LOCATION 
C  ***  IZS«Z  STARTING  LOCATION 
C  ***  IZE-Z  ENDING  LOCATION 
C  ***  IP-PLANE  X.  Y  OR  Z  <1.2  OR  3 > 

C  ***  IL-LOCATION  OF  PLANE 

C  «**#*#*##**»*»»»»»*»•*»***»»#**•*****#**** 

SUBROUTINE  WHERE 
COMMON  /C3/IXA.  IYB.  IL.  IV,  IP 
COMMON  /C4/IXS,  IXE,  IYS,  I  YE,  IZS,  IZE 
25  WRITE <5,  24) 

WRITE  <  5,  33)  IXS,  IXE,  IYS,  I  YE,  IZS,  IZE 
24  FORMAT < '  MIN/MAX  POINTS  ALONG  XYZ  AXES  TO  BE  VIEWED 
1  <  IE:  9,  9,  1,  11,  1,  A)  ') 

33  FORMAT< IX,  614) 

WRITEC5,  34) 

34  FORMAT <  1 X.  * ) 

READ  (  5,  *.  ERR*25,  END«*69)  IXS,  IXE,  IYS,  I  YE,  IZS,  IZE 

69  WRITE<5,  70)  IP 

70  FORMAT <  '  SPECIFY  PLANE  MX"  OR  "Y"  OR  "Z"  ',A1,  *) 
READ<5,  71,  END=1000)  IP 

71  FORMAT ( 1 A1 ) 

IFdP.  EG.  0)  GO  TO  1001 
1000  IF<IP. EG.  'X')  GO  TO  200 
IF(IP.  EG.  'Y')  GO  TO  400 
IFdP.  EG.  '!')  CO  TO  300 
WRITE <5, 3332) 

3332  FORMAT < '  SPECIFY  "X"  "Y"  OR  "Z"  ',*) 

READ<5,  71,  END=1001 )  IP 

200  CONTINUE 

201  WRITE<5,  72)  IP,  IL 

72  FORMAT (  *  LOCATION  ALONG  ',A1,  '-AXIS  ',12,  *> 

READ<  5,  1 99,  END«*74  )  IL 

199  FORMAT  < 1 14 ) 

IF<  ( IL.  LT.  1 ) .  OR.  <IL.  GT.  18)  )G0  TO  201 
IYB-JZE 
IXA-IYE 
CO  TO  74 

WRITE<5,  72)  IP,  IL 
READ <5,  199.  END«74)  IL 
IF<  <  IL.  OT.  6).  OR.  <  IL.  LT.  1>>  GO  TO  300 
IXA-IXE 
IYB- I YE 
GO  TO  74 


300 


o  o  o  n 
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400  WRITE(5,  72)  IP.  IL 

READ<5. 199, END=74)  IL 

IF(  (IL.  LT.  1).  OR.  (IL.  GT.  11 )  )  GO  TO  400 

IXA-IXE 

IYB=IZE 

74  WRITE(5, 26)  IV 

26  FORMAT (  '  VECTORS  XY,  XZ,  OR  YZ: 

1  TYPE  1,2  OR  3  RESPECTIVELY  Ml,  IX,  *> 

READ<  5.  *.  END=1001  )  IV 
1001  RETURN 
END 

C  ***************************************** 

C  THIS  SUBROUTINE  DOES  THE  ACTUAL  DRAWING  OF  THE 
C  VECTORS 

C  ****************************************** 

SUBROUTINE  PL0T2<IIX,  IIY,  I  DATA 1,  IDATA2 ) 

INTEGER*4  IMX.  IMY,  IDATA1 <18,  11 ),  IDATA2< 18,  11 ) 

COMMON  /C5/IGD,  ISCX,  ISCY,  ISKX,  ISKY,  IX,  IY 
DO  1  1=1, IIY 
IYY=IY+IGD*I 
DO  2  J=l,  I IX 
I XX=IX+IGD*J 

IMX*IDATA1 ( J,  I ) / ISCX  +  ISKX 
IMY=IDATA2( J,  I ) 7 ISCY+ISKY 
CALL  MOVABS ( I XX,  IYY ) 

CALL  DRWREL< IMX,  IMY) 

2  CONTINUE 

1  CONTINUE 

RETURN 
END 

C  **************************************** 

.THIS  SUBROUTINE  ASSIGNS  THE  APPROPRIATE  VALUES 
.  .  FROM  THE  THREE-DIMENSIONAL  ARRAY,  IDATA  TO  THE 
.  .  TWO-DIMENSIONAL  ARRAYS,  IDATA1  AND  IDATA2  FOR 
.  .  SUBSEQUENT  PLOTTING  USING  SUBROUTINE  PL0T2. 

C  ************************************************ 

SUBROUTINE  DATA( JJ,  I, J,  IDATA1,  IDATA2) 

INTEGER*4  I  DATA  ( 18,  6,  11 ),  IZ,  IZ1,  IZ2,  IZH 
INTECER*4  IDATA1 ( 18, 11), IDATA2(18,  11) 

COMMON  /C3/IXA,  IYB,  IL,  IV,  IP 

COMMON  /Cl/IHRD,  ITEST, IU,  INVELS, INVELE,  INVEL 
COMMON  /C4/IXS,  IXE,  IYS,  I  YE,  IZS,  IZE 

READ(  lU'JJ)  (  (  (IDATA  (  I,  J,  K),  I-IXS,  IXE),  IZS.  IZE),  K=  IYS,  I  YE) 
749  IF  ( IP.  EQ.  'X')  CO  TO  700 
IF(  IP.  EQ.  'Y' )  CO  TO  701 


n  n  o  ci  n  n  a  o  oo 


.  THE  FOLLOWING  DO  LOOPS  ARE  FOR  THE  —  Z  —  PLANE 

DO  750  J*l,  I YD 
DO  751  1-1. IXA 
IZ-IDATAd,  IL,  J) 

CALL  DECODE( IV.  IZ,  IZ1,  IZ2) 

IF<J.  EG.  1)  IZH-IZ2 
IDATA1 (I,  J)*IZl-499 
IDATA1 < I,  J)»- IDATA1 < I.  J) 

IDATA2< I* J)*IZ2-499 
IDATA2(  I.  J)  —  IDATA2<  I.  J) 

751  CONTINUE 

750  CONTINUE 

I-I-l 
J=J-1 
RETURN 

.  THE  FOLLOWING  DO  LOOPS  ARE  FOR  THE  —  X  —  PLANE 

700  DO  760  J*l, IXA 

DO  761  1*1, IYB 
I Z*IDATA < IL,  I,  J) 

CALL  DECODE <  I V,  IZ,  IZ2>  I Z 1 ) 

IF(U.  EQ.  1)  I ZH=  I  Z2 
IDATA1 < I,  J)»IZl-499 
IDATAKI,  J)*-IDATA1(I,  J) 

IDATA21I,  J)*IZ2-499 
IDATA2( I,  J)*-IDATA2< I,  J) 

761  CONTINUE 

760  CONTINUE 

1*1-1 
J.J-1 

RETURN 

.  THE  FOLLOWING  DO  LOOPS  ARE  FOR  THE  —  Y  —  PLANE 

DO  770  J-l# IYB 
DO  771  1*1, IXA 
I Z*IDATA< I,  J,  IL) 

CALL  DECODE< IV,  IZ.  IZ1, IZ2) 

IF < J.  EG.  DIZH-IZ2 
IDATAKI,  J)*IZl-499 
IDATAKI,  J)— IDATAKI,  U) 

IDATA2( I,  J)*IZ2-499 
IDATA2< I,  J)*-IDATA2( I.  J) 

CONTINUE 
CONTINUE 
1*1-1 
J*J-1 
RETURN 
END 


701 

C 

771 

770 


o  o 
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C  **************************************** 

C  THIS  SUBROUTINE  DECODES  THE  DATA  FROM  ITS 
C  STORED  9-DIGIT  INTEGER  FORM. 

C  ***************************************** 


SUBROUTINE  DECODEUV,  12,  121.  IZ2) 
INTEGERS  IZ. IZ1. IZ2, IZ3, IZHH 
IZ1*IZ/1E6 

IZ2«(IZ-IZ1*1E6>/1E3 
IF(  IV.  EQ.  1  )RETURN 
IZHH-IZ2 

I Z3=IZ-IZ1*1E6-I Z2*1E3 
I 12=123 

IF( IV  EQ. 2)  RETURN 

IZ1*IZHH 

RETURN 

END 


***************************************** 

C  THIS  SUBROUTINE  DRAWS  THE  SCHEMATIC  OF  THE 
C  TUNNEL. 

C  ****************************************** 

SUBROUTINE  TUNNEL 

COMMON  /C8/IXHP. IYHP, IHOL 

I XHP-680 

IYHP»490 

IHOL«l 

CALL  INIMATU.  1,  1,  50,  0.  0,  0) 

MAG*24 

DO  10  II»0.  90,  15 

I X*IFIX (MAG«COS< FLOAT < 1 1 > *3.  14159/180.  )> 
I Y*IFI X ( MAG*S IN <  FLOAT  < 1 1 ) *3.  14159/180.  )) 
CALL  M0VE3< IX,  IY,  0) 

CALL  DRAW3I-IX,  -IY.  0) 

CALL  M0VE3<  IX.  -IV,  0) 

CALL  DRAW3(-IX,  IY, 0) 

10  CONTINUE 
I X-200 
IY-140 
IZ-120 
NX— 180 
NZ— 60 
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CALL  M0VE3<+IX,  +  IY,  +  IZ) 
CALL  DRAW3<+IX, -IV, +  IZ> 
CALL  DR  AW3 <  NX , - 1 Y,  ♦  I Z  > 
CALL  DRAW3<NX,  +  IY.  +  IZ) 
CALL  DR AU3 <  + 1 X ,  + 1 Y ,  *11) 

CALL  DRAU3 <  +IX, ♦IV,  NZ) 
CALL  DRAW3<+IX, -IY,  NZ  ) 
CALL  DR AW3 ( ♦  I X »  - 1 Y ,  + 1 Z  > 

CALL  MOVE3<  +  IX, -IY. NZ) 
CALL  DRAW3<NX,  -I  Y»  NZ  ) 

CALL  DRAW3<NX.  +  IY.  NZ  > 

CALL  DRAW3<  +  IX, +IY, NZ) 

CALL  M0VE3<NX, +  IY, NZ ) 

CALL  DRAW3<NX. +IY, +IZ) 

CALL  M0VE3 ( NX , -IY. +IZ) 

CALL  DRAW3<NX, -IY,  NZ) 

IX=-160 

DO  1  1  =  1.  IB 

CALL  MOVE3< IX,  -IY,  IZ) 

CALL  DRAW3< IX, -IY,  (IZ-10)) 
CALL  DRAW3< IX, -IY,  < IZ  +  10) > 
CALL  M0VE3 ( I X ,  - 1 Y ,  ( I Z  +20 > ) 
CALL  ANMODE 
WR  ITE<  5,  20)1 
FORMAT < IX,  12) 

IX=IX+20 

CALL  M0VE3(0»  -IY,  < IZ+43) ) 
CALL  ANMODE 
WRITE<  5,  21 ) 

FORMAT ('  "X"  ') 


IX»200 

IY—100 

DO  2  1-1, 11 

CALL  M0VE3UX,  IY,  IZ) 

CALL  DRAW3< <IX-10), IY, IZ) 
CALL  DRAW3< < IX+1 0) , IY, IZ) 

CALL  M0VE3UIX+20),  <IY+8),  IZ) 
CALL  ANMODE 
WRITE( 5,  20)1 
IY-IY+20 

CALL  M0VE3( <IX+35),  0,  IZ) 

CALL  ANMODE 
WRITE (5, 22) 

FORMAT <'  "V"  ') 


r>  o 


IY— 140 
IX— 180 
12-80 

DO  5  I-l,  6.  1 

CALL  M0VE3UX,  I Y,  IZ) 

CALL  DRAW3<  < IX+10) , IY.  IZ) 

CALL  DRAW3<  (IX-10),  IY,  IZ) 

CALL  M0VE3< < IX+20),  IY,  (IZ-30)) 

CALL  ANMODE 
WRITE (5,  20)  I 
5  IZ-IZ-20 

CALL  M0VE3<  < IX-35) ,  IY,  0) 

CALL  ANMODE 
WR ITE( 5. 23) 

23  FORMAT (y  "Z"  ') 

RETURN 

END 

«*#*«-*********"»***«"»••*»******»***#*#***»**** 
SUBROUTINE  ORIGIN*  V) 

COMMON  /C5/ICD.  IS«,  »SY,  ISKX,  ISKY,  IX,  I 
COMMON  /Cl/IHKD, ITEST, IU, INVEL1, INVEL 
COMMON  /C3/IXA,  1 Y0,  IL,  IV,  IP 
IDEC-JJ*5 

IF  <  IP.  EQ.  '  Y ' ; CO  TO  10 
IF ( IP.  EQ  'Z')GO  TO  10 
DO  1  1*1, 11 
I YY*I Y+IGD*I 
IYY-IYY-4 
DO  2  J-1,6 
I XX-IX+IGD*J 
I XX-I XX-2 

CALL  MOVABS ( I XX,  IYY) 

CALL  ANCH0(43) 

2  CONTINUE 

1  CONTINUE 

GO  TO  20 

10  DO  11  1*1, IYB 
I YY-IY+TGD*I 
IYY-IYY-4 
DO  12  J-l, IXA 
I XX-I X+IGD*J 
I XX-I XX-2 

CALL  MOVABS ( I XX,  IYY) 

CALL  ANCH0<43) 

12  CONTINUE 

11  CONTINUE 


pj  < 


CONTINUE 

IF < .  NOT.  I HRD) RETURN 
CALL  M0VABS<770,  200) 

CALL  CZAXIS<0> 

CALL  ANMODE 
WRITE<  5»  3)  IP>  IL 

FORMAT (  '  PLOT  OF  A1 ,  '-PLANE  AT  LOCATION  ', 
CALL  M0VABS<770,  190) 

CALL  ANMODE 
WRITE<  5.  4)  I  DEO 

FORMAT <'  AT  BLADE  ANGLE  OF  '.13,  '  DEGREES') 

RETURN 

END 


r 


VOLTS. FTN 
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c  **#*#***#*«*#****#******#**###****#*####** 

C  THE  PROCRAM  VOLTS.  FTN  IS  A  DATA  AGUISITION 
C  PROGRAM  WHICH  CREATES  A  FILE  F0R006.  DAT 
C  CONTAINING  THE  VOLTAGES  FOR  THE  AVERAGE  OF 
C  THE  14  BLADE  REVOLUTIONS. 

C  THESE  ARE  AVA,  AVB,  AND  AVC  IN  THE  PROGRAM. 

C  THE  PROGRAM  OPERATES  ON  THE  DATA  STORED  IN 
C  THE  DIRECT  ACCESS  FILE  BUFFER.  DAT.  BUFFER 
C  .  DAT  IS  CREATED  USING  THE  PROGRAM  MAKEBUF. 

C  FTN.  THUS  ONCE  BUFFER  HAS  BEEN  CREATED 
C  THERE  IS  NO  NEED  TO  USE  THE  A/D  CONVERTER 
C  FOR  THIS  PROGRAM. 

C  ONCE  YOU  HAVE  F0R006.  DAT  USE  EASY  GRAPHING 
C  TO  DO  THE  PLOTTING.  YOU  MUST  FIRST  COPY 
C  FOR 006. DAT  TO  PT.  DAT,  THEN  TYPE  EZC  AND 
C  THEN  RUN  VEL.  OF  COURSE  YOU  WILL  NEED  TO 
C  CHANGE  THE  SCALES  AND  LABELS. 

C  ****##•»*•*■**■*****#**•***#*#*#*#**#****##*#**# 

INTEGER*4  ARRAY <  72 ) ,  NLOC, NVEL,  IVX,  IVY,  IVZ 

1  ,  I,  IX.  IY.  IZ 

INTEGER  BUFFER <1024,  3) 

C  ***  THE  ARRAYS  EAANG. EBANG,  AND  ECANG  ARE  STORING  THE 
C  ***  VELOCITIES  FOR  THE  14  REVS. 

DIMENSION  EAANG (14), EBANG  < 14 ) ,  ECANG ( 14 > ,  V(3),  U(3) 
1  ,  ASIN( 3) , AA ( 72 ) , SUM ( 14 ) ,  BB ( 72 ) ,  CC ( 72 ) 

1-1 
IND-1 

OPEN ( UNI T-l.  NAME* 'BUFFER.  DAT ', TYPE* 'OLD  ' , 
l  ACCESS-  'DIRECT') 

READU'l)  <BUFFER<N,  1  >,  N-l,  1024) 

READ< 1  '2>  (BUFFER(N, 2), N-l,  1024) 

READ( 1  '3)  (BUFFER (N, 3), N-l, 1024) 

WRITE (5, 5000) 

FORMAT <3X. '  BUFFER  WAS  FILLED') 


5000 


Ill 


DO  83  J-1.72 
C  ( 72-360/ 5) 

DO  9  M-1,3 
I A-J 

DO  8  JJ-1. 14 

C  14-REVOLUTIONS 

IF(M.  EQ.  1)  EAANO<JJ) -FLOAT (BUFFER (I A,  1  )  )/ (32768.  /8.  ) 
IF(M.  Ed.  2)  EBANC(JJ) -FLOAT (BUFFER (I A.  2> )/ (32768.  /8.  ) 
IF(M.  EQ.  3)  ECANC (JJ> -FLOAT ( BUFFER ( I A.  3)  )/ (32768.  /8.  ) 
IA-IA+72 

8  CONTINUE 

IF(M.  EQ.  1)  CALL  AVERAG( EAANG. AVA) 

IF(M.  EQ.  2)  CALL  AVER  AG  ( EB ANG,  AVB  > 

IF(M.  EG.  3)  CALL  AVER  AG  ( ECANC,  AVC) 

IF(M.  EQ.  1)  AA ( J )  —  AVA 
IF(M.  EQ.  2)  BB(J)-AVB 
IF(M.  EQ.  3)  CC(J)»AVC 
WRITE(5.  *) AVA# AA(J) 

9  CONTINUE 
C 

85  CONTINUE 

C 

990  I-I+l 

WRITE(  5.  489) 

489  FORMAT <4X. 'DO  YOU  WANT  AVA,  AVB,  OR  AVC?') 

READ (5. 567 ) ANSW 
567  FORMAT <1 A3) 

IF(ANSW.  EQ.  'AVB'  )  WRITE(6,*)BB 
IF(ANSW.  EQ.  'AVC')  WRITE(6,  *)CC 
WRITE(6,  *)AA 
WRITE (5,  *)  AA 
DO  1001  I U- 1 #  3 
1001  CLOSE (UNIT-IU) 

CLOSE (UN IT-6) 

STOP 

END 

SUBROUTINE  AVERAQ( A*  AV) 

DIMENSION  A( 14).  SUM( 14) 

SUM ( 1 ) -A  < 1 ) 

DO  9  1-1,  13 

9  SUM( 1+1 )-SUM( I )+A( 1+1 > 

A  V-SUM  ( 1 4 )  / 1 4. 

RETURN 

END 


112 


APPENDIX  D 

FILMED  ANIMATION  OF  VECTOR  PLOTS 

In  order  to  more  graphically  illustrate  the  instantaneous  nature 
of  the  data  of  Figures  4-4  through  4-6  it  was  decided  to  film  the 
vector  plots  as  they  appear  on  the  screen  of  a  TEKTRONIX  4014  graphics 
computer  terminal.  At  the  time  of  this  writing  seven  planes  of  data 
have  been  filmed;  the  X-Y  planes  at  Z/R  =0.2  and  0.4;  the  X-Z  planes 
at  Y/R  =  0.0,  -0.8  and  -0.53;  and  the  Y-Z  planes  at  X/R  =  -0.8  and 
-1.07.  Each  plane  shows  the  flow  every  5°  of  rotor  rotation  from  0° 
to  355°. 

The  dim  light  condition  made  it  necessary  to  film  each  "instan¬ 
taneous"  plot  in  the  single  frame  exposure  mode.  A  16  mm  camera  was 
used  with  a  shutter  speed  of  1/4  second  per  frame  focused  at  a  distance 
of  2  feet  from  the  screen.  Twenty- five  frames  were  exposed  for  each 
"instantaneous  plot"  which  when  projected  at  24  frames  per  second  gives 
a  viewing  time  of  about  1  second  per  plot.  Thus  one  rotor  revolution 
is  projected  in  about  75  seconds. 


